Downloaded from http://biorxiv.org/ on September 18, 2014 



bioRviv 

/Vbeta 

THE PREPRINT SERVER FOR BIOLOGY 



Response to climate change is related to species traits in the 
Lagomorpha 

Katie Leach, Ruth Kelly, Alison Cameron, et al. 
bioRxiv first posted online January 14, 2014 

Access the most recent version at doi: http://dx.doi.org/10.1101/001826 



Copyright The copyright holder for this preprint is the author/funder. All rights reserved. No 
reuse allowed without permission. 



Downloaded from http://biorxiv.org/ on September 18, 2014 



l Title: Response to climate change is related to species traits in the Lagomorpha 

2 

3 Authors: Katie Leach 1 *, Ruth Kelly 1 , Alison Cameron 2 ' 3 , W. Ian Montgomery 1 ' 2 and 

4 Neil Reid 1 

5 Affiliations: 

6 

7 l Quercus, School of Biological Sciences, Queen's University Belfast, Belfast, BT9 

8 7BL. Northern Ireland. 

9 School of Biological Sciences, Queen's University Belfast, Belfast, BT9 7BL. 

10 Northern Ireland. 

1 1 "Institute for Global Food Security (IGFS), Queen's University Belfast, Belfast, BT9 

12 5BN. Northern Ireland. 

13 

14 * Correspondence to: kleachOl @qub.ac.uk 
15 



Downloaded from http://biorxiv.org/ on September 18, 2014 

16 Abstract: Climate change during the last five decades has impacted on natural systems 

17 significantly and the rate of current climate change is of great concern among 

18 conservation biologists. Here, we assess the projected change in the bioclimatic 

19 envelopes of all 87 species in the mammalian order Lagomorpha under future climate 

20 using expertly validated species distribution models. Results suggest that climate 

21 change will impact more than two-thirds of Lagomorphs, with leporids (rabbits, hares 

22 and jackrabbits) likely to undertake poleward shifts with little overall change in range 

23 extent, whilst pikas are likely to show extreme shifts to higher altitudes associated with 

24 marked range declines, including the likely extinction of Kozlov's Pika (Ochotona 

25 koslowi). Species traits were associated with predictions of change, with smaller-bodied 

26 species more likely to exhibit range contractions and elevational increases, but showing 

27 little poleward movement. Lagomorphs vulnerable to climate change require urgent 

28 conservation management to mitigate range declines and/or extinctions. 
29 
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30 Main Text: The Lagomorpha is one of the most important mammalian orders 

31 economically and scientifically; as a major food resource, laboratory animals, valued 

32 game animals, pests of agricultural significance and the base of many predator/prey 

33 relationships providing insights into entire trophic systems (1). Species Distribution 

34 Models (SDMs) are widely used in ecology and relate species occurrences at known 

35 locations to environmental variables to produce models of environmental suitability, 

36 which can be spatially or temporally extrapolated to unexplored areas or into past or 

37 future conditions (2, 3). Lagomorphs are native or introduced on all continents except 

38 Antarctica, occur from sea level to >5 ,000m and from the equator to 80°N spanning a 

39 huge range of environmental conditions. Whilst other studies have modelled large 

40 numbers of species (4) or, for example, a few key species within each mammalian order 

41 (5), Lagomorphs due to their restricted diversity provide an opportunity to examine the 

42 response of every species yielding a detailed picture of change within an entire order for 

43 the first time. Furthermore, the small number of Lagomorph species, compared to other 

44 mammalian groups, means that datasets can be verified in detail, modelled individually 

45 and outputs expertly validated. 

46 We assessed projected change in the bioclimatic envelopes of all 87 lagomorph 

47 species under future climate change using a combination of climatic, topographical and 

48 environmental variables (<5). The taxonomy of the group is in constant flux but all 

49 belong to just two families: the Ochotonidae and the Leporidae. The Ochotonidae 

50 consists of one monotypic group in the genus Ochotona containing 25 species of small, 

51 social, high-altitude pikas. The Leporidae has 32 species of large, solitary, cursorial 

52 hares and jackrabbits in the genus Lepus and 30 species of medium-sized, semi-social, 

53 fossorial rabbits from 10 genera (1). A quarter of Lagomorphs are listed in the IUCN 
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54 Red List of Threatened Species (www.iucnredlist.org) with a notable number of highly 

55 range-restricted species, for example, the Tres Marias cottontail (Sylvilagus graysoni) 

1 2 

56 and black jackrabbit (Lepus insularis) which occupy islands <500km and ~80km 

57 respectively in the Gulf of California. In addition, pikas as high-altitude specialists are 

58 extremely susceptible to changes in their environment (7) with only small changes in 

59 temperature leading to heat stress (8). 

60 We used a framework for model validation based jointly on subjective expert 

61 evaluation and objective model evaluation statistics (Fig. SI). Expert validation, 

62 whereby an acknowledged expert on each species judges model predictions for current 

63 and past distributions (9), can be a useful tool prior to making future extrapolations. 

64 Forty- six Lagomorph experts, including 20 members of the Lagomorph Specialist 

65 Group (LSG) of the IUCN Species Survival Commission (SSC), were paired to one or 

66 more species and asked to assess whether model projections accurately, roughly or did 

67 not capture the current and past range of each species (good, medium and poor 

68 respectively) (<5). There was a high degree of agreement between expert classification 

69 and model Kappa (k) values (Fig. S2). k is an objective measure of prediction accuracy 

70 based on input species records and background points adjusted for the proportion of 

71 correct predictions expected by random chance (10). Fifty-eight species (67%) were 

72 deemed 'modellable' with expert validation classed as medium or good and k>0A and 

73 29 species (33%) rejected as 'unmodellable' with expert validation classed as poor 

74 and/or k<0A. Unmodellable species were 4 times more likely to be listed by the IUCN 

75 as Data Deficient than modellable species, with 8 unmodellable species (28%) listed as 

76 threatened. Hereafter, all results are for modellable species only and, therefore, are a 

77 conservative estimate of the impact of climate change on the order as a whole. 
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78 Changes in climate are predicted to have strong influences on the ecology and 

79 distribution of species (11, 12) with pronounced impacts on terrestrial biodiversity (4). 

80 Global changes in predicted Lagomorph species richness suggest that almost a third of 

81 the Earth's terrestrial surface (3 1 .5 million km ) is predicted to experience loss of 

82 Lagomorph species by 2080 (Fig. 1). The desert regions of North Eastern China and 

83 hills of Sichuan become increasingly unsuitable, losing up to 10 species by 2080, 

84 including the woolly hare (Lepus oiostolus) and Glover's pika (Ochotona gloveri); 

85 which are predicted to undergo dramatic movements to higher elevations. In contrast, 

86 species gains were notable across: (a) northern Eurasia, due to poleward movements of 

87 the mountain hare (Lepus timidus) whose boreo-alpine niche is highly sensitive to 

88 changes in temperature (13, 14) and (b) North America, where some regions e.g. the 

89 Upper Missouri catchment of Montana and North Dakota gain up to 5 species including 

90 the desert and eastern cottontails (Sylvilagus audubonii and S. floridanus respectively) 

91 which are predicted to exhibit strong poleward movements. The majority of Lagomorph 

92 species in Africa were classed as unmodellable and as such Figure 1 is data deficient for 

93 this continent. Parts of Africa are expected to become drier and warmer under future 

94 climate, with large increases in arid land (15) leading to negative consequences for 

95 Lagomorphs. Nevertheless, the results for modellable species suggest species gains in 

96 Southern Africa, due to a latitudinal shift in the range of the African savannah hare 

97 (Lepus microtis). 

98 Species are predicted to adapt their bioclimatic niche to future conditions, migrate to 

99 maintain their current niche, or become range restricted and undergo population decline, 

100 local extirpation or total extinction under future scenarios (16). Thirty-six Lagomorph 

101 species exhibited range loss (63%), 48 poleward movements (83%) and 51 elevational 
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102 increases (88%). Thirty-five species (60%) exhibited range declines and either poleward 

103 movements or elevational increases. On average, all three groups of Lagomorphs 

104 exhibited significant poleward shifts (^=1,228= 30.1, p<0.001) and elevational increases 

105 (F df=1 ,228= 19.7, p<0.001) by 2080 (Fig. 2, Table SI). The average poleward shift for the 

106 Order was 1.1° and elevational shift 165m; these values are much higher than averages 

107 calculated in a recent meta-analysis collating information on a wide variety of 

108 taxonomic groups (77). Pikas exhibited the most substantial mean increase in elevation 

109 becoming increasingly isolated on mountain summits (e.g. the Rockies in North 

1 10 America and the Tibetan Plateau and high Himalayas) resulting in a significant 31% 

1 1 1 range contraction (Fdf=i,228 =5.0, p=0.03). Pikas are known to be extremely susceptible 

112 to heat stress (#); thus, increases in temperature or unpredictable seasonality are likely 

113 to induce elevational rather than latitudinal shifts in these high-altitude specialists. 

114 The leporids, typically being lowland species, exhibited less substantial increases in 

1 15 elevation but greater poleward shifts. As global temperatures increase, northern latitudes 

116 will become more climatically suitable for southern leporids and, therefore, species 

117 bioclimatic envelopes will track poleward to match. Thus, rabbits, hares and jackrabbits 

118 all exhibited little overall change in the total extent of their ranges. The majority of the 

119 modellable rabbit species were from the genus Sylvilagus inhabiting south and eastern 

120 USA and Mexico, which are projected to become generally drier (18) inducing 

121 latitudinal shifts in species to track suitable habitats or vegetation. Thus, by shifting 

122 their ranges poleward, leporids are predicted to be able to maintain or increase their 

123 range extent suggesting they are considerably less sensitive to projected climate change 

124 than pikas. 
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125 Recent studies have found that species' morphological and life history traits can be 

126 good predictors of climate change responses (19). Phylogenetic generalised least 

127 squares (PGLS) regressions estimated lambda (X) to be close to zero (/?>0.05) for 

128 changes in range, poleward movement and elevation indicating that observed trends 

129 were independent of phylogeny (20). Generally, the members of each group were 

130 capable of showing a variety of responses i.e. species of pika, rabbits, hares or 

131 jackrabbits exhibited both increases and decreases in each response variable within each 

132 group (Fig. 2). There was a significant positive relationship between range change and 

133 adult body mass (/?=0.258, F d f = 4 !52 =2.308, p=0.021) with hares and jackrabbits generally 

134 increasing their range by 2080 and pikas exhibiting range contraction (Table S2; Fig. 

135 S3). This could occur due to a relationship between adult body mass and dispersal 

136 distance, which was used to clip model projections, but a Spearman's Rank correlation 

137 suggested no correlation between these variables in the dataset used. Adult body mass 

138 was positively associated with poleward movement (fi=0. 196, ^=5,49=1-989, p=0.047) 

139 and inversely related to elevational change (/?= -0. 183, ^=2,50=2.019, p=0.043). Smaller 

140 species, principally pikas and some rabbits, were typically more capable of elevational 

141 shifts in range whilst larger species, principally hares and jackrabbits, had a greater 

142 tendency for poleward shifts. This relationship is probably due to the large number of 

143 small-bodied species that dwell in mountainous regions (n=24) which have more 

144 opportunity to shift altitudinally. The average adult body mass for Lagomorphs living in 

145 mountainous regions was 836g, compared to 1.8kg in lowland areas. A lm elevational 

146 shift is equivalent to a 6km latitudinal shift (11), so it is easier for smaller species to 

147 shift altitudinally. The number of litters a species produces per year was positively 

148 related to latitudinal and maximum elevational shifts (/?=0.215, Fdf=5,49=2.731, p=0.006 
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149 and /?=0.160, Fdf = 2,53= 1.746, p=0.081 respectively) with more fecund species being 

150 capable of more extreme upward or poleward shifts, as in Moritz et al, (2008) (21). 

151 There was also a positive relationship between species dietary niche and the degree of 

152 poleward shift predicted (0=0.181, F df=5 , 4 9=2.190, p=0.029). 

153 The uncertainties in SDM using projected future climate scenarios are well 

154 described 6 . Models are vulnerable to sampling bias (22), spatial scale issues (25), lack 

155 of data for rare species (24), uncertainty regarding future climate conditions (25) and 

156 insufficient independent evaluation (23). We have tackled these explicitly by accounting 

157 for sampling bias by restricting the set of background points used, using data with the 

158 highest spatial resolution available (30 arc-second) and selecting species records to 

159 match, bootstrapping models for rare species with few records, averaging climatic data 

160 from five global circulation models and using a framework of joint expert and metric - 

161 driven model validation to segregate outputs into those that were modellable and 

162 unmodellable before subsequent analysis. 

163 Regardless of individual model outcomes (see Supplementary Materials), the overall 

164 trends observed across the order Lagomorpha as a whole are compelling. This study did 

165 not take account of shifts in habitats, vegetation or human impact in response to climate 

166 change, but we have shown that adaptation to future climate conditions may be possible 

167 as some species were predicted to exhibit poleward movements, with only modest shifts 

168 in range or elevation e.g. snowshoe hare (Lepus americanus) and Appalachian cottontail 

169 (Sylvilagus obscurus). 

170 The predicted changes in climate conditions are likely to have greater impacts on 

171 isolated Lagomorph populations, i.e. those on islands and in high-altitude systems. 

172 Lagomorph species occupying islands (n=6) will, on average, lose 8km of their ranges 
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2 

173 compared to 0.2km gain for continental species (n=52), whilst mountain dwelling 

2 2 

174 species (n=24) will lose 37km of their ranges compared to 25km" gain for lowland 

175 species (n=34). If changes continue at the rate currently predicted until the 2080s, then 

176 there may be no climatically suitable range available for some montane or isolated 

177 species e.g. the Tres Marias cottontail (Sylvilagus graysoni), black jackrabb it (Lepus 

178 insularis) or Gansu pika (Ochotona cansus). Conservation strategies, such as assisted 

179 migration, could be the only option for these highly range-restricted species. 

180 Furthermore, conservation management will need to focus on small-bodied mammals as 

181 these are predicted to show more dramatic responses to changing climate. Small 

182 mammals are key in food webs sustaining predator populations, impacting plant 

183 communities by grazing and soil biology and hydrology by burrowing. Thus, 

184 fundamental shifts in Lagomorphs globally will cause trophic cascade effects, especially 

185 in northern latitudes such as the cyclic systems of the Arctic. The advancing knowledge 

186 of past extinction rates for Lagomorphs (26), along with significantly better bioclimatic 

187 envelope modelling, could aid the prediction and prevention of future extinctions. Our 

188 models suggest that Kozlov's pika (Ochotona koslowi) may become extinct by the 

189 2080s as the elevational increases required to maintain its current bioclimatic envelope 

190 disappear as it reaches the maximum elevation available. Assessment of vulnerability to 

191 climate change is needed urgently to identify how and to what extent species, taxa, 

192 communities and ecosystems are susceptible to future changes, taking into account 

193 likely changes to vegetation, human pressures and interspecific interactions, and to 

194 direct conservation management in an efficient and effective manner. Multi-species 

195 approaches are likely to lead to more effective mitigation measures and contribute to 
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196 our understanding of the general principles underpinning the biogeographical and 

197 ecological consequences of climate change impacts. 
198 
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342 




360 Figure 1. Change in predicted Lagomorph species richness from the 1930s to 

361 2080s. A, Global patterns in predicted species loss and gain showing details in B, North 

362 America and C, Asia. White areas indicate no change in species richness or areas that 

363 are occupied by unmodellable species with uncertain outcomes (i.e. data deficient). 

364 
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406 Figure 2. Characterisation of predicted Lagomorph bioclimatic envelope change. 

407 Scatterplots show the linear relationship between range change (%) and A, poleward 

408 movement (°), B, elevational change (m) and C, poleward movement and elevational 

409 change. The numbers of species in each quadrant that exhibited positive or negative 

410 change on each axis are shown in plain text with percentages in bold parentheses. 

411 Temporal trends for D, range change, E, poleward movement and F, elevational change 
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412 + 1 standard error for each species group; pikas (red), rabbits (green) and hares and 

413 jackrabbits (black). 
414 
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423 

424 Materials and Methods: 

425 Species data 

426 A total of 139,686 records including all 87 Lagomorph species were either downloaded 

427 from the Global Biodiversity Information Facility (GBIF) Data Portal (data.gbif.org), 

428 collated from species experts or members of the IUCN Lagomorph Specialist Group 

429 and/or searching the literature. Taxonomic, spatial and temporal data accuracy and 

430 sampling bias issues were dealt with leaving 41,874 records of which 3,207 were pre- 

431 1950 and 38,667 were post- 1950. All records were checked against the latest IUCN 

432 taxonomy, if names did not match after studying taxonomic synonyms and previous 

433 names they were rejected, and if records for the target species fell outside the extent of 

434 the IUCN geographic range polygon they were also rejected. Records with spatial error 

435 estimates >2km were removed and duplicates were eliminated unless they occurred in 

436 different temporal periods (pre-1950 and post-1950). To remove spatial bias the 

437 background data was only selected from sites at which any lagomorph had been 

438 recorded, and to remove temporal bias the background data (10,000 points) had the 

439 same proportion of past to current points as there were past to current species records. 
440 
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441 Environmental parameterisation 

442 Empirical climate data were downloaded from WorldClim and the CCAFS GCM Data 

443 Portal at 30 arc-second resolution (~lkm grid cells). Records pre-1950 were associated 

444 with mean climate data from 1900-1949 and those post- 1950 were associated with mean 

445 data from 1950-2000. Projected future climate data were obtained from the IPCC 4 th 

446 Assessment Report using the A2 emissions scenario. Although the A2 future climate 

447 change scenario was originally described as 'extreme climate change' it now appears to 

448 best represent the trend in observed climate. Pierce et al. (2009) (27) report that using 

449 data averaged across five global circulation models (GCMs) is substantially better than 

450 any one individual model and significantly reduces model error. We used the following 

451 five GCMs: CCCma-CGCIvB. 1/T47, CNRM-CM3, CSIRO-Mk3.0, HadCIvB and 

452 NASA-GISS-ER. Future projections adopted the time periods: 2020s (2010-2039), 

453 2050s (2040-2069) and 2080s (2070-2099). 

454 Fifteen environmental variables were used in their raw format: mean annual 

455 temperature (°C), temperature seasonality (°C), mean temperature of the hottest month 

456 (°C), mean temperature of the coldest month (°C), mean annual precipitation (mm), 

457 mean precipitation of the wettest month (mm), mean precipitation of the driest month 

458 (mm) and precipitation seasonality (mm). Further composite variables were calculated 

459 including: annual evapotranspiration (mm), annual water balance (mm), the number of 

460 months with a positive water balance, mean annual Normalised Difference Vegetation 

461 Index (NDVI), human influence index, solar radiation and surface roughness index. 

462 Evapotranspiration was calculated using the Hargreaves equation (28): 
463 

464 Evapotranspiration = 0.0023 x (T mean + 17.8) x (T max - T m i n ) 0 5 x R a 
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465 where, 
466 

467 R a = ((24 x 60)/7t) x G sc x d r [(co s x sirup x sinS) + (coscp x cosS x sinooj]) x 0.408 

468 G sc = solar constant = 0.0820 MJ m-2 min- 1 

469 d r = inverse relative distance Earth-Sun = 1 + 0.033 cos ((27i/365) x J) 

470 co s = sunset hour angle [rad] = arcos [-tan (cp) tan (5)] 

471 (p = latitude [radians] (grid file in decimal degrees converted to radians (multiply by 

472 71/180)) 

473 5 = solar decimation [rad] = 0.409 sin (((2ti/365) x J) - 1 .39) 

474 J = number of days in the year. J at the middle of the month is approximately given by J 

475 = INTEGER (30.4 Month - 15) = 15 on average. 
476 

477 Annual evapotranspiration was taken as the sum of all monthly values. Annual water 

478 balance was calculated by subtracting annual evapotranspiration from mean annual 

479 precipitation. The number of months with a positive water balance was calculated by 

480 subtracting each monthly evapotranspiration from its corresponding monthly 

481 precipitation, and then converting these into a binary format, where a value greater than 

482 zero was given a value of one and a value less than zero was kept at zero (29). The 

483 twelve binary files were then summed to calculate the number of months with a positive 

484 water balance. 

485 Mean annual Normalised Difference Vegetation Index (NDVI) was calculated from 

486 monthly values which were downloaded from the EDIT Geoplatform 

487 (http://edit.csic.es/Soil-Vegetation-LandCover.html). NDVI is commonly used to 

488 measure the density of plant growth and is obtained from NOAA AVHRR satellite 
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489 images. A negative value indicates snow or ice, a value around 0 indicates barren areas, 

490 values between 0.2 and 0.3 indicate grassland, and values near 1 indicate rainforests 

491 (30). Human influence index was downloaded from the SEDAC website 

492 (http://sedac.ciesin.columbia.edu/data/set/wildareas-v2-human-influence-index- 

493 geographic). This was a composite of human population density, railways, roads, 

494 navigable rivers, coastlines, night-time lights, built-up areas and agricultural and urban 

495 land cover. Values within the index range from 0 to 64, where zero equalled no human 

496 influence and 64 represented maximum human influence (31). Solar radiation was 

497 calculated using the Spatial Analyst function in ArcGIS 10.1 (ESRI, California, USA). 

498 Solar radiation is defined as the total amount of incoming solar insolation (direct and 

499 diffuse), or global radiation, and was measured in watt hours per square meter or 

500 WH/m (32). An index of surface roughness was also calculated by finding the 

501 difference between maximum and minimum gradient values, based on a global Digital 

502 Elevation Model at 30 arc- second resolution downloaded from WorldClim (33). 

503 Altitude was not included as a variable independently because organisms perceive 

504 climatic and habitat variables as proxies for altitude (34). 
505 

506 Species distribution modelling 

507 SDMs were run using MaxEnt (version 3.3.3k) (35, 36). Models used two different sets 

508 of input data: i) pre- and post-1950s data, or, ii) post-1950s data only. For some species, 

509 the set of input data used affected the model outcome substantially. The 'species with 

510 data' (SWD) format was used for data entry pairing pre- 1950 records with mean 

511 climatic variables from 1900-1949 and post-1950 records with mean data from 1950- 

512 2000. The models were validated using either 10 replicate bootstrapping for species 
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513 with low numbers of records (<30) or 4-fold cross-validation for species with high 

514 numbers of records (>30). Linear, quadratic and product feature types were used. The 

515 10 percentile threshold was applied to define likely presence and absence of each 

516 species. Records for each species were associated with global land cover data 

517 downloaded from the ESA GlobCover 2009 Project and any land class not occupied by 

518 the target species was marked as unsuitable. Thus, two sets of model outputs were 

519 created: i) showing suitable bioclimatic space (i.e. model predictions), and, ii) 

520 bioclimatic space restricted to occupied habitat (i.e. model predictions clipped to 

521 suitable habitats). Outputs were then clipped to a minimum convex polygon of the 

522 species occurrence records, buffered by a distance value specific to each species which 

523 was verified by experts from the IUCN Lagomorph Specialist Group, to remove areas 

524 of over-prediction. This meant that the suitable area within reach of each species was 

525 predicted, and led to much more realistic outputs than just simply predicting suitable 

526 area. 

527 Introduced ranges for the European hare (Lepus europaeus) and rabbit {Oryctolagus 

528 cuniculus) were not modelled due to lack of records (common species are frequently 

529 under-recorded). Mountain hare {Lepus timidus) populations in Ireland and mainland 

530 Eurasia were modelled separately due to the distinct ecological differences between the 

531 two (37) but the outputs were subsequently combined to reflect the current classification 

532 of the Irish hare (L. t. hibernicus) as a sub-species of L. timidus. 
533 

534 Model evaluation 

535 A bespoke website (lagomorphclimatechange.wordpress.com) was created to allow each 

536 expert to review the output of their allocated species. Experts selected either: i) pre- and 
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537 post-1950 input data showing the suitable bioclimatic envelope, ii) pre- and post-1950 

538 input data showing the suitable bioclimatic envelope restricted to suitable habitat, iii) 

539 post- 1950 input data only showing the suitable bioclimatic envelope, or, iv) post- 1950 

540 input data only showing the suitable bioclimatic envelope restricted to suitable habitat. 

541 Experts selected the model perceived as optimal in describing the known pre- 1950 and 

542 post- 1950 range of species overlaid by the current IUCN range extent polygon. Expert 

543 evaluation classed each model as 'good', 'medium' or 'poor' according to the criteria in 

544 Anderson et al. (2003) (9). In addition to expert evaluation, models were assessed using 

545 Kappa as calculated in R (version 3.0.0). The 'accuracy' function in the SDMTools 

546 package (38) was used to calculate the Area under the Curve (AUC) value, omission 

547 rate, sensitivity, specificity, proportion correct, Kappa and True Skill Statistic (TSS). 

548 Models that had a Kappa >0.4 and those that were ranked as either 'good' or 'medium' 

549 by expert evaluation were defined as modellable and were carried forward for analysis, 

550 whilst those that had a Kappa <0.4 or were ranked as 'poor' by expert evaluation were 

551 defined as unmodellable and were rejected from further analysis (see Fig. SI). Kappa 

552 was used because it utilises commission and omission errors (39), and although it does 

553 not take into account prevalence like TSS (40), it has widely accepted thresholds (41, 

554 42, 43) which make it best suited for use in this study. The optimum threshold for 

555 Kappa was taken as 0.4 (40, 41, 42). The model settings selected as providing the 

556 optimal outcome were used to predict species bioclimatic envelopes under future 

557 climate during the 2020s, 2050s and 2080s. Future predictions were clipped to the 

558 buffered minimum convex polygon of the target species further buffered by the 

559 dispersal distance (kilometres/year) of each species multiplied by the number of years 

560 elapsed from the present (1950-2000) taken as 1975. Annual dispersal distances were 
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561 elicited from each species expert during the model evaluation procedure or taken from 

562 the literature. For some species where no data were available, the average for specific 

563 groups (i.e. Asian pikas or African hares) was used. Predicted range size, mean latitude 

564 and minimum, mean and maximum elevation for each species and each time period 

565 were calculated. Model outputs for each species can be found in the Supplementary 

566 Materials. 
567 

568 Species traits 

569 Species trait data were downloaded from the PanTHERIA database (44) and updated by 

570 searching the literature. Correlated traits were removed to reduce collinearity and the 

571 final set of traits used to describe each modellable species were: activity cycle 

572 (nocturnal only, diurnal only, mixed), adult body mass (g), diet breadth (number of 

573 dietary categories eaten from: vertebrates, invertebrates, fruit, flowers/nectar/pollen, 

574 leaves/branches/bark, seeds, grass and roots/tubers), gestation length (days), habitat 

575 breadth (above ground-dwelling, aquatic, fossorial and ground-dwelling), home range 

2 2 

576 (km ), litters per year, litter size, population density (n/km ) and age at sexual maturity 

577 (days). 
578 

579 Statistical analysis 

580 To illustrate global change in the distribution of Lagomorphs, the difference in 

581 predicted species richness per cell was calculated between the 1930s (1900-1949) and 

582 2080s (2070-2099). The difference in model output metrics (range size, mean latitude 

583 and minimum, mean and maximum elevation) was calculated between the 1930s and 

584 2080s. Change in range size was expressed in percentage change but change in latitude 
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585 was represented as degree movements towards poles and change in elevation in metres. 

586 GLMs were used to test the differences in temporal trends for range change, poleward 

587 movements and elevation between lagomorph groups: pikas, rabbits and hares and 

588 jackrabbits. PGLS regressions were then performed using the 'caper' (45) package in R 

589 to test whether changes in model predictions varied with species traits. A lagomorph 

590 phylogeny was extracted from the mammalian supertree provided by Fritz et al. (2009) 

591 (46). Likely clade membership for five species not included in this phylogeny was 

592 determined from Ge et al. (2013) (26) and the literature, and then missing tips were 

593 grafted on using an expanded tree approach (47). Outliers were removed prior to 

594 analysis identified as those with a studentized residual >3 following Cooper et al. 

595 (2012) (48). All subsequent models exhibited normally distributed residuals tested using 

596 Shapiro-Wilk. The scaling parameter lambda (k) varies from 0, where traits are 

597 independent of phylogeny, to 1, where species' traits covary in direct proportion to their 

598 shared evolutionary history (20). We estimated X for the overall model and tested its 

599 significance during the PGLS analysis. All subset regressions were run using the 

600 'dredge' function in the 'MuMIn' (49) package in R and model averaging was used to 

601 describe the effect of each variate. 
602 
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YES 



Expert evaluation 
= Medium / Good 
■ 



Kappa value 
>0.4 



NO 



Expert evaluation 
= Poor 
■ 



YES 



NO 
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NO 



MODELLABLE 



I- UNMODEL LAB L E 



UNMODELLABLE 



UNMODEL LAB L E 



Figure SI. Framework for assessing whether species were modellable or 
unmodellable based on Kappa values and expert evaluation classification. The 

optimum threshold for Kappa was taken as 0.4 (40, 41, 42). Expert evaluations were 
classified according to Anderson et al (2003) (9). 
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611 

iModellable (retained) Unmodellable (rejected) 




Good Medium Poor 



612 Expert Evaluation 

613 Figure S2. Agreement between expert evaluation and model accuracy. Mean Kappa 

614 +1 standard error within the categories assigned by expert evaluation. Dark grey bars 

615 indicate species that were deemed modellable and retained for further analyses, whereas 

616 light grey bras indicate unmodellable species that were rejected. Sample sizes (i.e. 

617 numbers of species) are shown in the bars. 
618 
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Diet breadth Number of litters per year Adult body mass (g) 



Figure S3. Relationships between species traits and responses to future climate 
change. The ability of species' traits to predict changes in A, range B, mean elevation 
C, maximum elevation and D, E, F, poleward movement under future climate (between 
-1930s and ~2080s) for each group; pikas (red), rabbits (green) and hares and 
jackrabbits (black). Diet breadth is a categorical variable and is therefore represented as 
a bar plot (+1 standard error) with sample sizes (i.e. numbers of species) shown above 
the bars. Only significant predictors of change are shown here. The dashed line at zero 
indicates no change in the response variable. 
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630 Table SI. Results of GLMs characterising predicted Lagomorph bioclimatic 

631 envelope changes. Significant p values are in bold. Group refers to Lagomorph 

632 taxonomy, i.e. pikas, rabbits and hares & jackrabbits. 

633 



Response variable 


Term 


Deviance 


Residual 


F 


Df 


P 








Deviance 








Range change (km) 


Group 


0.217 


4.334 


6.467 


2,229 


0.002 




Year 


0.085 


4.249 


5.033 


1,228 


0.026 




Group: Year 


0.055 


4.195 


1.639 


2,226 


0.196 


Mean latitudinal change 


Group 


<0.001 


<0.001 


0.000 


2,229 


1.000 


0 


Year 


<0.001 


<0.001 


30.102 


1,228 


<0.001 




Group: Year 


<0.001 


<0.001 


8.375 


2,226 


<0.001 


Mean elevation change 


Group 


1442316 


9159926 


19.856 


2,9 


<0.001 


(m) 


Year 


716859 


8443066 


19.737 


1, 8 


<0.001 




Group: Year 


234776 


8208290 


3.232 


2,6 


0.041 



634 
635 
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636 Table S2. Results of phylogenetically controlled generalised least square 

637 regressions. Significant p values for model- averaged coefficients are in bold. F and p 

638 values for the top model are listed under each response variable; asterisks (*) indicate 

639 traits in the top model. Lambda (X) confidence intervals and significance from 0 and 1 

640 are also shown. 
641 



Response variable 


Trait 


P ± s.e. 


F 


P 


Range change (km) 


Adult body mass (g)* 


0.258 ±0.112 


2.308 


0.021 


F if=4 , 52=4.28, />=0.005 


Diet breadth 


0.137 ±0.088 


1.552 


0.121 




Gestation length (days) 


-0.154 ±0.121 


1.269 


0.204 


Lambda 95%CI=0, 0.147 


Litters per year 


0.090 ± 0.084 


1.073 


0.283 


N.S. from 0,p<0.01 from 1 


Litter size 


-0.062 + 0.098 


0.634 


0.526 




Home range size (km 2 ) 


-0.069 + 0.099 


0.698 


0.485 




Population density (n/km 2 ) 


-0.079 ± 0.090 


0.878 


0.380 




Age at sexual maturity 


-0.062 ± 0.089 


0.694 


0.488 




Activity cycle 


0.020 ±0.134 


0.151 


0.880 




Habitat breadth 


0.050 ±0.101 


0.493 


0.622 


Mean latitudinal change (°) 


Adult body mass (g)* 


0.196 ±0.099 


1.989 


0.047 


^df=5,49=6.10, p<0.001 


Diet breadth* 


0.181 ±0.082 


2.190 


0.029 




Litter size* 


0.128 ±0.097 


1.320 


0.187 


Lambda 95%CI=0, 0.209 


Litters per year* 


0.215 ±0.079 


2.731 


0.006 


N.S. from 0,p<0.01 from 1 


Activity cycle 


0.097 ±0.124 


0.787 


0.431 




Age at sexual maturity 


-0.106 ±0.084 


1.254 


0.210 




Home range size (km 2 ) 


-0.088 ±0.113 


0.787 


0.438 




Gestation length (days) 


0.147 ±0.132 


1.112 


0.266 




Habitat breadth 


0.028 + 0.093 


0.304 


0.761 




Population density (n/km 2 ) 


0.027 ± 0.084 


0.318 


0.750 


Mean elevation change (m) 


Adult body mass (g)* 


-0.183 ± 0.091 


2.019 


0.043 


F if=2 , 50=5.92, p=0.005 


Gestation length (days) 


0.099 ±0.106 


0.932 


0.351 




Diet breadth 


-0.1 10 ±0.079 


1.386 


0.166 


Lambda 95%CI=0, 0.205 


Home range size (km 2 ) 


0.074 ± 0.075 


0.986 


0.324 


N.S. from 0,p<0.01 from 1 


Litters per year 


-0.057 ± 0.070 


0.813 


0.416 




Age at sexual maturity 


0.055 ± 0.075 


0.731 


0.324 




Activity cycle 


0.075 ±0.105 


0.708 


0.479 




Population density (n/km 2 ) 


0.064 + 0.071 


0.902 


0.367 




Litter size 


0.019 ±0.079 


0.235 


0.814 




Habitat breadth 


-0.005 ± 0.085 


0.058 


0.953 


Max. elevation change (m) 


Litters per year * 


0.160 ±0.092 


1.746 


0.081 


F df=2 , 53=3.54, p=0.036 


Gestation length (days) 


-0.136 ±0.116 


1.169 


0.242 




Habitat breadth 


0.095 ±0.102 


0.934 


0.351 


Lambda 95%CI=0, 0.384 


Age at sexual maturity 


-0.098 + 0.096 


1.020 


0.308 


N.S. from 0,p<0.01 from 1 


Activity cycle 


0.041 ±0.118 


0.348 


0.728 




Litter size 


0.014 ±0.104 


0.137 


0.891 




Population density (n/km 2 ) 


0.054 ±0.091 


0.596 


0.551 




Adult body mass (g) 


0.034 ±0.121 


0.280 


0.780 




Home range size (km 2 ) 


-0.003 ±0.107 


0.029 


0.977 




Diet breadth 


0.004 ± 0.093 


0.047 


0.962 



642 
643 
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644 Species accounts (#1-89) 

645 

646 Key for response curves 



Abbreviation 


Description 


HII 


Human Influence index 


Hilliness 


Surface Roughness index 


MaxPrec 


Maximum precipitation 


MaxTemp 


Maximum temperature 


MinPrec 


Minimum precipitation 


MinTemp 


Minimum temperature 


PrecSea 


Precipitation seasonality 


RealMAR 


Mean annual precipitation 


RealMAT 


Mean annual temperature 


TempSea 


Temperature seasonality 


etpsum 


Annual evapotranspiration 


ndvi 


Normalised difference vegetation index 


radiation 


Solar radiation 


wbann 


Annual water balance 


wbpos 


Number of months with a positive water balance 



647 
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Model evaluation metric 



Species account #1 - Pygmy rabbit {Brachylagus idahoensis) 
n = 39 

Expert: Penny Becker, Washington Dept. of Fish & Wildlife, USA 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 15km/year (Expert) 

Status: MODELLABLE 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


1.00 


Proportion correct 


1.00 


Kappa 


0.75 


True Skill Statistic 


0.90 



Summary: The Pygmy rabbit's bioclimatic envelope is predicted to decline by 87% with a 1° 
mean latitudinal poleward shift and mean increase in elevation of ~300m driven predominately 
by an increase in mean minimum elevation (>600m) with little change in mean maximum 
elevation (~50m). 95% of the permutation importance of the model was contributed to by mean 
annual temperature (64.5%), maximum temperature (25.2%) and annual water balance (5.9%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Beauvais GP, et al. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#2 - Riverine rabbit (Bunolagus monticularis) 
n= 109 

Status: Kai Collins, University of Pretoria, South Africa 

Expert evaluation: Good 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 7.5km/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.98 


Omission rate 


0.03 


Sensitivity 


0.97 


Specificity 


1.00 


Proportion correct 


1.00 


Kappa 


0.85 


True Skill Statistic 


0.97 



Summary: The Riverine rabbit's bioclimatic envelope is predicted to decline by 85% with a 
-1° mean ^tjtiudinal poleward ^jft and mean increa^ in elevation of g-^0m driven by sjjnilar 
increases in both minimum and maximum elevation. 95% of the permutation importance of the 
model was contributed to by minimum temperature (33.1%) and precipitation (22.5%), 
temperature seasonality (22.3%), mean annual temperature (6.7%) and precipitation (5.0%), 
annual evapotranspiration (4.3%) and precipitation seasonality (2.6%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Beauvais GP, et al. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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Model evaluation metric 



#3 - Hispid hare (Caprolagus hispidus) 
n = 18 

Expert: Gopinathan Maheswaran, Zoological Survey of 
India 

Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: 5km/year (Expert) 
Status: MODELLABLE 



AUC 


0.97 


Omission rate 


0.06 


Sensitivity 


0.94 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.81 


True Skill Statistic 


0.94 



Summary: The Hispid hare's bioclimatic envelope is predicted to increase by 21% with a 
-1.5° mean latitudinal poleward shift and mean increase in elevation of ~70m driven by 
increases in maximum elevation. 95% of the permutation importance of the model was 
contributed to by mean annual temperature (52.7%), precipitation seasonality (29.0%), annual 
evapotranspiration (6.6%), number of months with a positive water balance (2.9%), maximum 
precipitation (2.2%) and minimum temperature (1.7%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Maheswaran, G. & Smith, A.T. 
2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#4 - Antelope jackrabbit (Lepus alleni) 
n = 32 

Expert: Paul Krausman, University of Montana 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: 25km/year (Expert) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.91 


Omission rate 


0.16 


Sensitivity 


0.84 


Specificity 


0.99 


Proportion correct 


0.98 


Kappa 


0.26 


True Skill Statistic 


0.83 



Summary: The Antelope jackrabbit' s bioclimatic envelope is predicted to increase by 172% 
with a -3° mean latitudinal poleward shift and mean increase in elevation of ~20m driven by 
increases in maximum elevation. 95% of the permutation importance of the model was 
contributed to by precipitation seasonality (24.8%), minimum precipitation (21.1%), annual 
water balance (16.3%), minimum temperature (13.0%), temperature seasonality (7.3%), 
normalised difference vegetation index (5.8%), annual evapotranspiration (5.2%) and human 
influence index (3.7%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Malpica, F.J. & Rangel Cordero, 
H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#5 - Snowshoe hare (Lepus americanus) 
n = 506 

Expert: Charles Krebs, University of British Colombia & 

Rudy Boonstra, University of Toronto 

Expert evaluation: Good 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: 24km/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.07 


Sensitivity 


0.93 


Specificity 


0.97 


Proportion correct 


0.97 


Kappa 


0.72 


True Skill Statistic 


0.90 



Summary: The Snowshoe hare's bioclimatic envelope is predicted to decline by 7% with a 
~2°mean latitudinal poleward shift and mean decrease in elevation of ~10m, but with increases 
in both minimum and maximum elevation. 95% of the permutation importance of the model 
was contributed to by mean annual temperature (83.1%), maximum temperature (7.2%), 
normalised difference vegetation index (2.1%), mean annual precipitation (1.5%) and annual 
evapotranspiration (1.2%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Murray, D. & Smith, A.T. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#6 - Arctic hare (Lepus arcticus) 
n = 18 

Expert: David Gray, Grayhound Information Services 
Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: 2km/year (Chapman & Flux, 1990) 
Status: UNMODELLABLE 

Summary: The Arctic hare's bioclimatic envelope is predicted to decline by 30% with a -2° 
mean latitudinal poleward shift and mean decrease in elevation of ~10m driven by decreases in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
maximum temperature (34.0%), annual evapotranspiration (15.4%), surface roughness index 
(14.4%), human influence index (10.1%), normalised difference vegetation index (9.2%), mean 
annual temperature (4.4%), precipitation seasonality (1.5%), mean annual precipitation (1.5%), 
maximum precipitation (0.3%), minimum precipitation (0.3%), solar radiation (0.3%), 
minimum temperature (0. 1 %) and number of months with a positive water balance (0. 1 %). 




a) -1930s (1900-1949) b) ~ 1980s (1950-2000) c) ~ 2020s (2010-2039) 



Model evaluation metric 



AUC 


0.97 


Omission rate 


0.06 


Sensitivity 


0.94 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.36 


True Skill Statistic 


0.94 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Murray, D. & Smith, AT. 
2008) 
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#7 - Japanese hare (Lepus brachyurus) 
n = 9 

Expert: Koji Shimano, Shinshu University, Japan 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: lkm/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.43 


True Skill Statistic 


0.99 



Summary: The Japanese hare's bioclimatic envelope is predicted to increase by 9% with no 
latitudinal poleward shift and a mean increase in elevation of ~10m driven by a decrease in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
temperature seasonality (31.8%), annual water balance (20.6%), human influence index 
(13.5%), mean annual precipitation (8.9%), maximum precipitation (6.2%), precipitation 
seasonality (5.4%) and minimum precipitation (4.4%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Yamada, F. & Smith, A.T. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#8 - Black-tailed jackrabbit (Lepus californicus) 

n = 970 

Expert: Alejandro Velasquez, UNAM-Canada 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 1.5km/year (Average for NA leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.93 


Omission rate 


0.07 


Sensitivity 


0.93 


Specificity 


0.94 


Proportion correct 


0.94 


Kappa 


0.69 


True Skill Statistic 


0.87 



Summary: The Black-tailed jackrabbit' s bioclimatic envelope is predicted to decline by 25% 
with a -2° mean latitudinal poleward shift and mean decrease in elevation of ~75m, but with 
increases in both minimum and maximum elevation. 95% of the permutation importance of the 
model was contributed to by precipitation seasonality (31.8%), annual evapotranspiration 
(18.5%), maximum temperature (17.0%), mean annual temperature (8.9%), minimum 
temperature (8.0%), minimum precipitation (4.0%), human influence index (3.9%) and 
temperature seasonality (3.5%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 



Black = predicted suitable 
envelope 




Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#9 - White-sided jackrabbit (Lepus callotis) 
n = 31 

Expert: Jennifer Frey, New Mexico State University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 18.9km/year (Average for NA leporids) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.97 


Omission rate 


0.05 


Sensitivity 


0.95 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.36 


True Skill Statistic 


0.93 



Summary: The White-sided jackrabbit' s bioclimatic envelope is predicted to increase by 3% 
with a ~l°mean latitudinal poleward shift and a mean increase in elevation of ~150m driven by 
an increase in maximum elevation. 95% of the permutation importance of the model was 
contributed to by precipitation seasonality (35.4%), annual evapotranspiration (22.3%), 
minimum temperature (17.0%), mean annual temperature (10.5%), minimum precipitation 
(6.5%), maximum precipitation (6.5%), surface roughness index (3.0%) and number of months 
with a positive water balance (2.7%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#10 - Cape hare (Lepus capensis) 
n = 231 

Expert: John Flux, IUCN Lagomorph Specialist Group 

Expert evaluation: Poor 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 35km/year (Expert) 

Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.93 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.97 


Proportion correct 


0.97 


Kappa 


0.56 


True Skill Statistic 


0.87 



Summary: The Cape hare's bioclimatic envelope is predicted to decrease by 45% with -2° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~330m driven by 
an increase in maximum elevation. 95% of the permutation importance of the model was 
contributed to by annual evapotranspiration (33.1%), minimum precipitation (29.6%), 
maximum temperature (9.7%), human influence index (7.2%), normalised difference vegetation 
index (4.6%), minimum temperature (3.2%), number of months with a positive water balance 
(2.9%), maximum precipitation (2.2%), mean annual precipitation (2.1%) and precipitation 
seasonality (2.0%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Drew, C, et al. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#11 - Broom hare (Lepus castroviejoi) 
n = 164 

Expert: Pelayo Acevedo, University of Porto 

Expert evaluation: Medium 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: lkm/year (Average for European 

leporids) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.11 


Sensitivity 


0.89 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.80 


True Skill Statistic 


0.89 



Summary: The Broom hare's bioclimatic envelope is predicted to decrease by 90% with a 
-0.2° mean latitudinal poleward shift and a mean increase in elevation of ~450m driven by an 
increase in minimum elevation. 95% of the permutation importance of the model was 
contributed to by mean annual temperature (62.0%), maximum temperature (20.6%), 
temperature seasonality (10.9%) and surface roughness index (3.8%). 




a) -1930s (1900-1949) b) ~ 1980s (1950-2000) c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#12 - Yunnan hare (Lepus comas) 
n = 59 

Expert: Weihe Yang, Institute of Zoology, Chinese 
Academy of Sciences 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 2.5km/year (Average for Asian leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.92 


Omission rate 


0.15 


Sensitivity 


0.85 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.67 


True Skill Statistic 


0.84 



Summary: The Yunnan hare's bioclimatic envelope is predicted to decrease by 65% with a -0.1° 
mean latitudinal poleward shift and a mean increase in elevation of ~100m driven by both increases 
in maximum and minimum elevation. 95% of the permutation importance of the model was 
contributed to by precipitation seasonality (59.9%), maximum temperature (26.8%), temperature 
seasonality (2.9%), number of months with a positive water balance (2.6%), annual 
evapotranspiration (2.0%) and mean annual temperature (1.7%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 
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#13 - Korean hare (Lepus coreanus) 
n = 6 

Expert: Weihe Yang, Institute of Zoology, Chinese 
Academy of Sciences 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 2.5km/year (Average for Asian leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.86 


True Skill Statistic 


0.99 



Summary: The Korean hare's bioclimatic envelope is predicted to increase by 500% with a -8° 
mean latitudinal poleward shift and a mean increase in elevation of ~70m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
temperature seasonality (27.2%), mean annual precipitation (25.6%), minimum temperature 
(17.7%), annual water balance (13.1%), normalised difference vegetation index (4.7%), 
precipitation seasonality (4.5%) and human influence index (2.4%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




d) ~ 2050s (2040-2069) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



2080s (2070-2099) 
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#14 - Apennine hare (Lepus corsicanus) 
n = 59 

Expert: Francesco Angelici, Italian Foundation of 

Vertebrate Zoology 

Expert evaluation: Medium 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.53 


True Skill Statistic 


0.99 



Summary: The Apennine hare's bioclimatic envelope is predicted to increase by 125% with a -2° 
mean latitudinal poleward shift and a mean decrease in elevation of ~60m. 95% of the permutation 
importance of the model was contributed to by minimum temperature (37.9%), annual 
evapotranspiration (34.3%), temperature seasonality (11.4%), minimum precipitation (11.1%) and 
maximum temperature (2.8%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Angelici, F.M., Randi, E., Riga, 
F. & Trocchi, V. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#15 - European hare (Lepus europaeus)-native range only 
n = 6,186 

Expert: Neil Reid, Queen's University Belfast 

Expert evaluation: Medium 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: 2km/year (Chapman & Flux, 1990) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.81 


Omission rate 


0.07 


Sensitivity 


0.93 


Specificity 


0.69 


Proportion correct 


0.78 


Kappa 


0.57 


True Skill Statistic 


0.62 



Summary: The European hare's bioclimatic envelope is predicted to increase by 50% with a -3° 
mean latitudinal poleward shift and a mean increase in elevation of ~20m driven by an increase in 
both maximum and minimum elevation. 95% of the permutation importance of the model was 
contributed to by annual evapotranspiration (58.0%), minimum temperature (9.0%), temperature 
seasonality (7.3%), surface roughness index (5.5%), human influence index (5.1%), precipitation 
seasonality (3.3%), annual water balance (3.3%), minimum precipitation (3.0%), normalised 
difference vegetation index (2.4%) and mean annual temperature (1.5%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 
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#16 - Ethiopian hare (Lepus fagani) 
n = 9 

Expert: Zelalem Tolesa, Addis Ababa University 
Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 8km/year (Average for African leporids) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.56 


True Skill Statistic 


0.99 



Summary: The Ethiopian hare's bioclimatic envelope is predicted to decrease by 15% with no 
latitudinal poleward shift and a mean increase in elevation of ~200m driven by an increase in 
maximum and minimum elevation. 95% of the permutation importance of the model was 
contributed to by temperature seasonality (92.7%) and annual evapotranspiration (5.3%). 




a) ~ 1930s (1900-1949) 




d) ~ 2050s (2040-2069) 




b) ~ 1980s (1950-2000) 




e) ~ 2080s (2070-2099) 



W0 r< 

c) ~ 2020s (2010-2039) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 




#17 - Tehuantepec jackrabbit (Lepus flavigularis) 
n = 8 

Expert: Arturo Carillo-Reyes, Universidad de Ciencias y 

Artes de Chiapas 

Expert evaluation: Poor 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Expert) 

Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.95 


True Skill Statistic 


0.99 



Summary: The Tehuantepec jackrabbit's bioclimatic envelope is predicted to decrease by 45% 
with a -1° mean latitudinal poleward shift and a mean increase in elevation of ~450m driven by an 
increase in maximum and minimum elevation. 95% of the permutation importance of the model 
was contributed to by temperature seasonality (81.7%), mean annual temperature (2.7%) and 
normalised difference vegetation index (1.7%). 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Cervantes, F.A., et al. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#18 - Iberian hare (Lepus granatensis) 
n = 1675 

Expert: Pelayo Acevedo, University of Porto 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: 7km/year (Expert) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.08 


Sensitivity 


0.92 


Specificity 


0.95 


Proportion correct 


0.95 


Kappa 


0.81 


True Skill Statistic 


0.87 



Summary: The Iberian hare's bioclimatic envelope is predicted to increase by 40% with a -1° 
mean latitudinal poleward shift and a mean increase in elevation of ~10m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
maximum precipitation (39.0%), annual evapotranspiration (38.0%), minimum temperature (15.0%) 
and maximum temperature (3.0%). 




a) -1930s (1900-1949) b) ~ 1980s (1950-2000) c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#19 - Abyssinian hare (Lepus habessinicus) 
n=l 

Expert: Zelalem Tolesa, Addis Ababa University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 25km/year (Average for African leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.82 


True Skill Statistic 


0.99 



Summary: The Abyssinian hare's bioclimatic envelope is predicted to decrease by 4% with a -4° 
mean latitudinal shift towards the Equator and a mean decrease in elevation of ~200m driven by an 
decrease in maximum elevation. 95% of the permutation importance of the model was contributed 
to by temperature seasonality (87.6%), annual evapotranspiration (4.9%) and minimum 
precipitation (4.1%). 





#20 - Hainan hare (Lepus hainanus) 
n = 9 

Expert: Youhua Chen, Wuhan University, China 
Expert evaluation: Good 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Average for island leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.86 


True Skill Statistic 


0.99 



Summary: The Hainan hare's bioclimatic envelope is predicted to increase by 4% with no 
latitudinal poleward shift and a mean increase in elevation of ~20m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
minimum temperature (71.0%), mean annual temperature (22.8%) and temperature seasonality 
(5.4%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Lazell, ].,etal. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#21 - Black jackrabbit (Lepus insularis) 

n = 3 

Expert: Tamara Rioja Pardela, Universidad de Ciencias y 
Artes de Chiapas, Mexico 
Expert evaluation: Good 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Average for island leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


1.00 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


1.00 


Proportion correct 


1.00 


Kappa 


1.00 


True Skill Statistic 


1.00 



Summary: The Black jackrabbit's bioclimatic envelope is predicted to decrease by 100% with a 
-0.3° mean latitudinal polewards shift and a mean increase in elevation of ~50m driven by an 
increase in both minimum and maximum elevation. 95% of the permutation importance of the 
model was contributed to by precipitation seasonality (28.6%), minimum precipitation (16.5%), 
annual water balance (13.3%), minimum temperature (12.9%), mean annual temperature (7.1%), 
temperature seasonality (5.5%), surface roughness index (5.0%) and normalised difference 
vegetation index (3.2%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 




#22 - Manchurian hare (Lepus mandshuricus) 
n = 36 

Expert: Deyan Ge, Institute of Zoology, Chinese Academy of 
Sciences 

Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Sokolov, V.E. et ah, 2009) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.96 


Omission rate 


0.08 


Sensitivity 


0.92 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.78 


True Skill Statistic 


0.92 



Summary: The Manchurian hare's bioclimatic envelope is predicted to decrease by 50% with a 
~l°mean latitudinal polewards shift and a mean increase in elevation of ~70m driven by an increase 
in maximum and minimum elevation. 95% of the permutation importance of the model was 
contributed to by precipitation seasonality (77.3%), mean annual temperature (10.4%), minimum 
temperature (4.5%), normalised difference vegetation index (2.6%) and annual evapotranspiration 
(2.4%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 





#23 - African savannah hare (Lepus microtis) 
n = 82 

Expert: John Flux, IUCN Lagomorph Specialist Group 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic only 

Dispersal distance: 15km/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.93 


Omission rate 


0.13 


Sensitivity 


0.87 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.62 


True Skill Statistic 


0.86 



Summary: The African savannah hare's bioclimatic envelope is predicted to decrease by 15% 
with a -3° mean latitudinal polewards shift and a mean increase in elevation of ~100m driven by an 
increase in maximum and minimum elevation. 95% of the permutation importance of the model 
was contributed to by mean annual temperature (40.9%), maximum temperature (18.5%), annual 
evapotranspiration (9.2%), minimum precipitation (6.2%), temperature seasonality (5.4%), mean 
annual precipitation (5.4%), minimum temperature (4.9%), precipitation seasonality (4.1%) and 
annual water balance (1.7%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#24 - Indian hare (Lepus nigricollis) 
n = ll 

Expert: Gopinathan Maheswaran, Zoological Survey of 
India 

Expert evaluation: Good 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: 6km/year (Expert) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.59 


True Skill Statistic 


0.99 



Summary: The Indian hare's bioclimatic envelope is predicted to decrease by 10% with a -2° 
mean latitudinal polewards shift and a mean increase in elevation of ~80m driven by an increase in 
maximum and minimum elevation. 95% of the permutation importance of the model was 
contributed to by precipitation seasonality (48.0%), mean annual temperature (14.9%), human 
influence index (11.9%), minimum temperature (10.3%), temperature seasonality (7.0%), number 
of months with a positive water balance (2.8%) and minimum precipitation (1.5%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Maheswaran, G. & Jordan, M. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#25 - Woolly hare (Lepus oiostolus) 
n = 84 

Expert: Weihe Yang, Institute of Zoology, Chinese 
Academy of Sciences 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 2.5km/year (Average for Asian leporids) 
Status: MODELLABLE 

Summary: The Woolly hare's bioclimatic envelope is predicted to decrease by 25% with a -1° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~680m driven by an 
increase in maximum and minimum elevation. 95% of the permutation importance of the model 
was contributed to by minimum precipitation (82.1%), maximum temperature (5.0%), minimum 
temperature (4.6%) and annual evapotranspiration (3.8%). 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.11 


Sensitivity 


0.89 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.63 


True Skill Statistic 


0.89 




a) ~ 1930s (1900-1949) b) ~ 1980s (1950-2000) c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Indian CAMP Workshop & 



Johnston, C.H. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#26 - Alaskan hare (Lepus othus) 
n = 8 

Expert: Eric Waltari, City University of New York 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic only 

Dispersal distance: 2km/year (Average for Arctic leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.89 


True Skill Statistic 


0.99 



Summary: The Alaskan hare's bioclimatic envelope is predicted to increase by 80% with a -3° 
mean latitudinal polewards shift and a mean increase in elevation of ~100m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by annual 
evapotranspiration (52.8%), mean annual temperature (26.4%), precipitation seasonality (11.3%) 
and mean annual precipitation (5.0%). 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Murray, D. & Smith, A.T. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#27 - Burmese hare (Lepus peguensis) 
n=l 

Expert: Thomas Gray, WWF Greater Mekong 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 2.5km/year (Average for Asian leporids) 
Status: MODELLABLE 



Model evaluation metric 

AUC 0.99 

Omission rate 0.00 

Sensitivity 1.00 

Specificity 0.99 

Proportion correct 0.99 

Kappa 0.58 

True Skill Statistic 0.99 



Summary: The Burmese hare's bioclimatic envelope is predicted to increase by 40% with a -2° 
mean latitudinal polewards shift and a mean increase in elevation of ~ 1 80m driven by an increase in 
minimum and maximum elevation. 95% of the permutation importance of the model was 
contributed to by temperature seasonality (52.6%), minimum precipitation (32.9%), solar radiation 
(6.5%) and mean annual temperature (4.3%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Duckworth, J.W., et oil. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




-D.389 D.983 25.526 3993013 -66.741 4334.332 



#28 - Scrub hare {Lepus saxatilis) 
n = 39 

Expert: Kai Collins, University of Pretoria 
Expert evaluation: Poor 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 25km/year (Average for African leporids) 
Status: UNMODELLABLE 



Model evaluation metric 

AUC 0.95 

Omission rate 0.08 

Sensitivity 0.92 

Specificity 0.97 

Proportion correct 0.97 

Kappa 0.18 

True Skill Statistic 0.89 



Summary: The Scrub hare's bioclimatic envelope is predicted to decrease by 15% with a -2° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~65m driven by an 
increase in maximum and minimum elevation. 95% of the permutation importance of the model 
was contributed to by annual evapotranspiration (52.8%), mean annual temperature (26.4%), 
precipitation seasonality (11.3%) and mean annual precipitation (5.0%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Collins, K., et al, 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#29 - Chinese hare (Lepus sinensis) 
n = 141 

Expert: Weihe Yang, Institute of Zoology, Chinese 
Academy of Sciences 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 2.5km/year (Average for Asian leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.11 


Sensitivity 


0.89 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.81 


True Skill Statistic 


0.89 



Summary: The Chinese hare's bioclimatic envelope is predicted to increase by 60% with a -2° 
mean latitudinal polewards shift and a mean increase in elevation of ~25m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by mean 
annual temperature (56.5%), temperature seasonality (22.9%), precipitation seasonality (8.0%), 
mean annual precipitation (5.4%) and annual evapotranspiration (3.5%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#30 - Ethiopian highland hare (Lepus starcki) 

n = 13 

Expert: Zelalem Tolesa, Addis Ababa University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 8km/year (Average for African leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.96 


Omission rate 


0.08 


Sensitivity 


0.92 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.86 


True Skill Statistic 


0.92 



Summary: The Ethiopian highland hare's bioclimatic envelope is predicted to decrease by 90% 
with a -7° mean latitudinal shift towards the Equator and a mean decrease in elevation of ~140m 
driven by a decrease in maximum elevation. 95% of the permutation importance of the model was 
contributed to by temperature seasonality (80.6%) and minimum temperature (18.5%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 




#31 - Desert hare (Lepus tibetanus) 
n = 55 

Expert: Chelmala Srinivasulu, Osmania University, India 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 2.5km/year (Average for Asian leporids) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.92 


Omission rate 


0.15 


Sensitivity 


0.85 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.57 


True Skill Statistic 


0.85 



Summary: The Desert hare's bioclimatic envelope is predicted to decrease by 50% with no 
latitudinal shift towards the Equator, but a mean increase in elevation of ~320m driven by an 
increase in maximum elevation. 95% of the permutation importance of the model was contributed 
to by minimum precipitation (54.5%), minimum temperature (19.3%), maximum precipitation 
(17.3%), annual evapotranspiration (2.8%) and human influence index (1.7%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(China Red List & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#32 - Mountain hare (Lepus timidus) - Eurasian ^ Model evaluation metric 

populations 
n = 2,460 

Expert: Neil Reid, Queen's University Belfast 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 
Dispersal distance: 2km/year (Expert) 
Status: MODELLABLE 



AUC 


0.91 


Omission rate 


0.08 


Sensitivity 


0.92 


Specificity 


0.90 


Proportion correct 


0.91 


Kappa 


0.74 


True Skill Statistic 


0.82 



Summary: The Mountain hare's bioclimatic envelope is predicted to decrease by 10% with a -4° 
mean latitudinal polewards shift and a mean decrease in elevation of ~10m driven by a decrease in 
maximum elevation. 95% of the permutation importance of the model was contributed to by annual 
evapotranspiration (87.9%), temperature seasonality (3.9%), minimum precipitation (2.1%) and 
minimum temperature (1.6%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 

d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#33 - Irish hare (Lepus timidus hibernicus) ^Model evaluation metric ^ 

n = 706 Omission rate 0.08 

Expert: Neil Reid, Queen's University Belfast Sensitivity 0.92 

Expert evaluation: Medium Specificity 0.97 

Data: Only modern Proportion correct 0.96 

Envelope: Climatic and habitat Treadll Statistic 0^88 
Dispersal distance: 2km/year (Expert) 
Status: MODELLABLE 



Summary: The Irish hare's bioclimatic envelope is predicted to decrease by 50% with a -0.5° 
mean latitudinal polewards shift and a mean decrease in elevation of ~10m driven by a decrease in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
temperature seasonality (44.6%), annual evapotranspiration (41.5%), normalised difference 
vegetation index (6.4%), precipitation seasonality (3.6%) and maximum precipitation (2.5%). 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#34 - Mountain hare (Lepus timidus) ^Model evaluation metric ^ 

- Eurasian & Irish populations combined Omission rate 0.01 

n= 3,166 Sensitivity 0.93 

Expert: Neil Reid, Queen's University Belfast Specificity 0.91 

Expert evaluation: Medium " Proportion correct 0.91 

Data: Only modern xSkill Statistic 0^84 
Envelope: Climatic and habitat 
Dispersal distance: 2km/year (Expert) 
Status: MODELLABLE 



Summary: The Mountain hare's bioclimatic envelope is predicted to decrease by 10% with a -2° 
mean latitudinal polewards shift and a mean decrease in elevation of ~40m driven by a decrease in 
maximum elevation. 95% of the permutation importance of the model was contributed to by annual 
evapotranspiration (87.6%), temperature seasonality (4.1%), minimum precipitation (2.1%) and 
minimum temperature (1.6%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 





#35 - Tolai hare (Lepus tolai) 
n = 316 

Expert: Chelmala Srinivasulu, Osmania University, India 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 2.5km/year (Average for Asian leporids) 
Status: MODELLABLE 

Summary: The Tolai hare's bioclimatic envelope is predicted to increase by 70% with a ~3°mean 
latitudinal polewards shift and a mean increase in elevation of ~280m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
temperature seasonality (43.5%), maximum temperature (17.3%), mean annual temperature 
(15.2%), precipitation seasonality (9.6%), annual evapotranspiration (4.6%), mean annual 
precipitation (3.7%) and minimum precipitation (1.8%). 




a) -1930s (1900-1949) b) ~ 1980s (1950-2000) c) ~ 2020s (2010-2039) 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.11 


Sensitivity 


0.89 


Specificity 


0.99 


Proportion correct 


0.98 


Kappa 


0.76 


True Skill Statistic 


0.88 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(China Red List & Johnston, C.H. 
2008) 
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#36 - White-tailed jackrabbit (Lepus townsendii) 

n = 275 

Expert: Eric Waltari, City University of New York 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 18.9km/year (Average for NA leporids) 
Status: MODELLABLE 

Summary: The White-tailed jackrabbit' s bioclimatic envelope is predicted to decrease by 10% 
with a ~4° mean latitudinal polewards shift and a mean increase in elevation of ~200m driven by an 
increase in minimum elevation. 95% of the permutation importance of the model was contributed to 
by mean annual temperature (62.5%), maximum temperature (28.5%), temperature seasonality 
(3.4%) and annual evapotranspiration (1.6%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 



Model evaluation metric 



AUC 0.94 

Omission rate 0.10 

Sensitivity 0.90 

Specificity 0.99 

Proportion correct 0.99 

Kappa 0.76 

True Skill Statistic 0.89 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, AT. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 





#37 - Yarkand hare {Lepus yarkandensis) 
n=49 

Expert: Weihe Yang, Institute of Zoology, Chinese 

Academy of Sciences 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 2km/year (Smith & Xie, 2008) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.74 


True Skill Statistic 


0.90 



Summary: The Yarkand hare's bioclimatic envelope is predicted to decrease by 100% with a -1° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~ 1000m driven by an 
increase in minimum elevation. 95% of the permutation importance of the model was contributed to 
by minimum precipitation (58.8%), maximum precipitation (16.3%), mean annual precipitation 
(11.9%) and minimum temperature (11.7%). 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#38 - Sumatran striped rabbit (Nesolagus netscheri) 
n = 11 

Expert: Hariyo Wibisono, Wildlife Conservation Society, 
Indonesia 

Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: O.Olkm/year (Expert) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.95 


True Skill Statistic 


0.99 



Summary: The Sumatran striped rabbit's bioclimatic envelope is predicted to decrease by 91% 
with a -1° mean latitudinal shift towards the Equator and a mean increase in elevation of ~330m 
driven by an increase in minimum elevation. 95% of the permutation importance of the model was 
contributed to by temperature seasonality (99.3%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Meijaard, E. & Sugardjito, J. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#39 - Annamite striped rabbit (Nesolagus timminsi) 
n=4 

Expert: Thomas Gray, WWF Greater Mekong & Andrew 

Tilker, University of Texas Austin 

Expert evaluation: Poor 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: lOkm/year (Expert) 

Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.50 


True Skill Statistic 


0.99 



Summary: The Annamite striped rabbit's bioclimatic envelope is predicted to decrease by 1500% 
with a -0.3° mean latitudinal shift towards the Equator and a mean decrease in elevation of ~20m 
driven by an decrease in minimum elevation. 95% of the permutation importance of the model was 
contributed to by solar radiation (43.3%), minimum temperature (30.0%), mean annual temperature 
(14.1%), human influence index (6.4%) and annual evapotranspiration (2.3%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Abramov, A., et al. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#40 - Alpine pika (Ochotona alpina) 
n = 16 

Expert: Sumiya Ganzorig, Hokkaido University 
Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: lOkm/year (Similar ecology to O.pallasi) 
Status: UNMODELLABLE 

Summary: The Alpine pika's bioclimatic envelope is predicted to decrease by 10% with a -4° 
mean latitudinal polewards shift and a mean decrease in elevation of ~80m. 95% of the permutation 
importance of the model was contributed to by minimum precipitation (50.1%), mean annual 
temperature (31.1%), number of months with a positive water balance (10.9%), minimum 
temperature (2.3%) and annual evapotranspiration (1.3%). 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.29 


True Skill Statistic 


0.99 






#41 - Silver pika (Ochotona argentata) 
n=4 

Expert: Andrew Smith, Arizona State University 

Expert evaluation: Poor 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Average for Asian pikas) 

Status: UNMODELLABLE 



Model evaluation metric 



AUC 


1.00 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


1.00 


Proportion correct 


1.00 


Kappa 


1.00 


True Skill Statistic 


1.00 



Summary: The Silver pika's bioclimatic envelope is predicted to decrease by 80% with a -4° 
mean latitudinal polewards shift and a mean increase in elevation of ~360m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
minimum precipitation (62.7%), human influence index (17.1%), mean annual precipitation (8.6%), 
minimum temperature (5.4%) and annual evapotranspiration (1.8%). 
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a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 






#42 - Gansu pika (Ochotona cansus) 
n = 38 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 1.5km/year (Similar ecology to O.roylei) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.11 


Sensitivity 


0.89 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.61 


True Skill Statistic 


0.89 



Summary: The Gansu pika's bioclimatic envelope is predicted to decrease by 60% with a -0.4° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~230m driven by an 
increase in minimum elevation. 95% of the permutation importance of the model was contributed to 
by minimum precipitation (89.8%) and minimum temperature (6.9%). 




a) ~ 1930s (1900-1949) 




d) ~ 2050s (2040-2069) 




b) ~ 1980s (1950-2000) 




e) ~ 2080s (2070-2099) 




c) ~ 2020s (2010-2039) 

Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 




#43 - Collared pika (Ochotona collaris) 
n = 193 

Expert: Hayley Lanier, University of Michigan & David 

Hik, University of Alberta 

Expert evaluation: Poor 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: lkm/year (Expert) 

Status: UNMODELLABLE 



Model evaluation metric 

AUC 0.95 

Omission rate 0.10 

Sensitivity 0.90 

Specificity 0.99 

Proportion correct 0.99 

Kappa 0.86 

True Skill Statistic 0.90 



Summary: The Collared pika's bioclimatic envelope is predicted to increase by 20% with a -2° 
mean latitudinal polewards shift and a mean decrease in elevation of ~140m. 95% of the 
permutation importance of the model was contributed to by annual evapotranspiration (86.7%), 
mean annual temperature (3.3%), normalised difference vegetation index (3.2%) and maximum 
temperature (2.0%). 
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a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 




#44 - Plateau pika (Ochotona curzoniae) 
n = 131 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Good 
Data: Only modern 
Envelope: Climatic and habitat 
Dispersal distance: O.lkm/year (Expert) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.11 


Sensitivity 


0.89 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.76 


True Skill Statistic 


0.88 



Summary: The Plateau pika's bioclimatic envelope is predicted to decrease by 30% with a -1° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~700m driven by an 
increase in minimum and maximum elevation. 95% of the permutation importance of the model 
was contributed to by minimum precipitation (48.0%), maximum temperature (42.0%) and annual 
evapotranspiration (6.8%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#45 - Daurian pika (Ochotona dauurica) 
n = 131 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: O.lkm/year (Similar ecology to O.curzoniae) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.66 


True Skill Statistic 


0.89 



Summary: The Daurian pika's bioclimatic envelope is predicted to decrease by 25% with a -1° 
mean latitudinal polewards shift and a mean decrease in elevation of ~60m. 95% of the permutation 
importance of the model was contributed to by minimum precipitation (92.3%) and number of 
months with a positive water balance (2.7%). 




a) -1930s (1900-1949) 




b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#46 - Chinese red pika (Ochotona erythrotis) 
n = 39 

Expert: Andrew Smith, Arizona State University 

Expert evaluation: Poor 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Average for Asian pikas) 

Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.53 


True Skill Statistic 


0.89 



Summary: The Chinese red pika's bioclimatic envelope is predicted to decrease by 20% with a 
-3° mean latitudinal polewards shift and a mean decrease in elevation of ~2400m driven by a 
decrease in minimum elevation. 95% of the permutation importance of the model was contributed 
to by minimum precipitation (87.2%), minimum temperature (7.2%) and mean annual temperature 
(2.9%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#47 - Forrest's pika (Ochotona forresti) 
n = 9 

Expert: Andrew Smith, Arizona State University 

Expert evaluation: Poor 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Average for Asian pikas) 

Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.53 


True Skill Statistic 


0.89 



Summary: The Forrest's pika's bioclimatic envelope is predicted to decrease by 40% with a -1° 
mean latitudinal polewards shift and a mean increase in elevation of ~600m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
minimum precipitation (36.8%), temperature seasonality (20.6%), normalised difference vegetation 
index (12.1%), maximum temperature (9.3%), precipitation seasonality (7.6%), human influence 
index (6.2%) and surface roughness index (4.2%). 








i 





a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#48 - Glover's pika (Ochotona gloveri) 
n = 22 

Expert: Andrew Smith, Arizona State University 

Expert evaluation: Medium 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Average for Asian pikas) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.65 


True Skill Statistic 


0.99 



Summary: The Glover's pika's bioclimatic envelope is predicted to decrease by 50% with a -0.5° 
mean latitudinal polewards shift and a mean increase in elevation of ~270m driven by an increase in 
minimum and maximum elevation. 95% of the permutation importance of the model was 
contributed to by minimum precipitation (46.3%), minimum temperature (28.8%), mean annual 
temperature (12.7%), human influence index (3.8%), temperature seasonality (3.1%) and 
precipitation seasonality (1.8%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#49 - Hoffmann's pika (Ochotona hoffmanni) 
n = 5 

Expert: Andrey Lissovsky, Zoological Museum of Moscow 

State University 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Average for Asian pikas) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.91 


True Skill Statistic 


0.99 



Summary: The Hoffmann's pika's bioclimatic envelope is predicted to decrease by 90% with a 
~l°mean latitudinal shift towards the Equator and a mean increase in elevation of ~470m driven by 
an increase in minimum elevation. 95% of the permutation importance of the model was 
contributed to by precipitation seasonality (73.2%), mean annual temperature (12.6%), minimum 
temperature (6.5%) and minimum precipitation (5.3%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#50 - Siberian pika {Ochotona hyperborea) 
n = 16 

Expert: Julia Witczuk, Warsaw Agricultural University, 
Poland 

Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: lOkm/year (Similar ecology to O.alpina) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.97 


Omission rate 


0.06 


Sensitivity 


0.94 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.40 


True Skill Statistic 


0.93 



Summary: The Siberian pika's bioclimatic envelope is predicted to decrease by 100% with a -5° 
mean latitudinal polewards shift and a mean increase in elevation of ~200m driven by an increase in 
minimum and maximum elevation. 95% of the permutation importance of the model was 
contributed to by mean annual temperature (72.5%), annual evapotranspiration (10.4%), minimum 
precipitation (6.4%), human influence index (2.8%), normalised difference vegetation index (2.6%) 
and annual water balance (1.8%). 



















1/ * 




a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#51- Hi pika (Ochotona iliensis) 
n = 11 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Poor 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 1 km/year (Similar ecology to O.koslowi) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.55 


True Skill Statistic 


0.99 



Summary: The Hi pika's bioclimatic envelope is predicted to decrease by 20% with a -1° mean 
latitudinal polewards shift and a mean increase in elevation of ~80m. 95% of the permutation 
importance of the model was contributed to by minimum precipitation (52.7%), minimum 
temperature (22.2%), mean annual temperature (8.8%), maximum precipitation (6.5%), temperature 
seasonality (3.0%) and human influence index (2.2%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#52 - Kozlov's pika {Ochotona koslowi) ^Model evaluation metric ^ 

n = 5 Omission rate 0.00 

Expert: Andrew Smith, Arizona State University Sensitivity 1.00 

Expert evaluation: Medium Specificity 0.99 

Data: Only modern Proportion correct 0.99 

Envelope: Climatic and habitat TrTskill Statistic 0^99 
Dispersal distance: lkm/year (Expert) 
Status: MODELLABLE 



Summary: The Kozlov's pika's bioclimatic envelope is predicted to decrease by 100% (total 
extinction). 95% of the permutation importance of the model was contributed to by mean annual 
precipitation (66.8%), minimum precipitation (24.2%) and human influence index (5.8%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 




#53 - Ladak pika (Ochotona ladacensis) 
n = 18 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 0.05km/year (Similar ecology to O.curzoniae) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.63 


True Skill Statistic 


0.99 



Summary: The Ladak pika's bioclimatic envelope is predicted to decrease by 70% with a -1° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~ 1400m driven by an 
increase in minimum and maximum elevation. 95% of the permutation importance of the model 
was contributed to by mean annual precipitation (47.9%), minimum precipitation (41.7%) and 
minimum temperature (5.7%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#54 - Large-eared pika (Ochotona macrotis) 
n =49 

Expert: Nishma Dahal, National Centre for Biological 

Sciences, India 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 1 km/year (Similar ecology to O.roylei) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.43 


True Skill Statistic 


0.89 



Summary: The Large-eared pika's bioclimatic envelope is predicted to decrease by 40% with a 
-1° mean latitudinal shift towards the Equator and a mean increase in elevation of ~880m driven by 
an increase in minimum elevation. 95% of the permutation importance of the model was 
contributed to by minimum precipitation (87.3%), minimum temperature (5.8%) and mean annual 
temperature (2.1%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#55 - Nubra's pika (Ochotona nubrica) 
n = 13 

Expert: Nishma Dahal, National Centre for Biological 

Sciences, India 

Expert evaluation: Medium 

Data: Only modern 

Envelope: Climatic and habitat 

Dispersal distance: 0.05km/year (Similar ecology to O.curzoniae) 
Status: UNMODELLABLE 



Model evaluation metric 

AUC 0.99 

Omission rate 0.00 

Sensitivity 1.00 

Specificity 0.99 

Proportion correct 0.99 

Kappa 0.35 

True Skill Statistic 0.99 



Summary: The Nubra's pika's bioclimatic envelope is predicted to increase by 1% with no 
latitudinal polewards shift, but a mean increase in elevation of ~200m driven by an increase in 
minimum and maximum elevation. 95% of the permutation importance of the model was 
contributed to by minimum precipitation (76.1%) and mean annual temperature (20.9%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A. F. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#56 - Pallas's pika (Ochotona pallasi) 
n = 19 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: lOkm/year (Sokolov, V.E. et al., 2009) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.46 


True Skill Statistic 


0.99 



Summary: The Pallas's pika's bioclimatic envelope is predicted to decrease by 60% with a -2° 
mean latitudinal polewards shift and a mean increase in elevation of ~40m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
minimum precipitation (46.1%), mean annual precipitation (35.5%), mean annual temperature 
(9.4%) and minimum temperature (6.7%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#57 - American pika (Ochotona princeps) 
n = 670 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 
Dispersal distance: 16.1km/year (Expert) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.9 


Specificity 


0.99 


Proportion correct 


0.98 


Kappa 


0.87 


True Skill Statistic 


0.89 



Summary: The American pika's bioclimatic envelope is predicted to decrease by 25% with a -1° 
mean latitudinal polewards shift and a mean decrease in elevation of ~10m driven by an decrease in 
minimum elevation. 95% of the permutation importance of the model was contributed to by mean 
annual temperature (88.1%), annual evapotranspiration (5.0%) and maximum temperature (2.8%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Beever, E. & Smith, A.T. 2011) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#58 - Little pika (Ochotona pusilla) 
n = 30 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 4km/year (Sokolov, V.E. et ah, 2009) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.93 


Omission rate 


0.13 


Sensitivity 


0.87 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.58 


True Skill Statistic 


0.86 



Summary: The Little pika's bioclimatic envelope is predicted to increase by 20% with a -1° mean 
latitudinal polewards shift and a mean increase in elevation of ~70m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
temperature seasonality (70.2%), annual water balance (13.3%), mean annual temperature (10.2%) 
and minimum temperature (4.5%). 




a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 
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#59 - Royle's pika (Ochotona roylei) 
n = 22 

Expert: Sabuj Bhattacharya, Wildlife Institute of India 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: lkm/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.98 


Omission rate 


0.05 


Sensitivity 


0.95 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.76 


True Skill Statistic 


0.95 



Summary: The Royle's pika's bioclimatic envelope is predicted to decrease by 20% with no 
latitudinal polewards shift and a mean increase in elevation of ~340m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
minimum precipitation (53.2%), maximum temperature (20.4%), temperature seasonality (10.2%), 
human influence index (5.5%), surface roughness index (4.2%) and precipitation seasonality 
(3.2%). 





#60 - Afghan pika (Ochotona rufescens) 
n = ll 

Expert: Chelmala Srinivasulu, Osmania University, India 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Average for Asian pikas) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.74 


True Skill Statistic 


0.99 



Summary: The Afghan pika's bioclimatic envelope is predicted to increase by 5% with a -1° 
mean latitudinal polewards shift and a mean increase in elevation of ~380m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
minimum precipitation (89.2%), minimum temperature (5.5%) and normalised difference 
vegetation index (2.3%). 




a) ~ 1930s (1900-1949) 



b)~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A.F. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#61 - Turkestan red pika (Ochotona rutila) 

n = 13 

Expert: Andrey Lissovsky, Zoological Museum of Moscow 

State University 

Expert evaluation: Poor 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Average for Asian pikas) 

Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.25 


True Skill Statistic 


0.99 



Summary: The Turkestan red pika's bioclimatic envelope is predicted to decrease by 10% with a 
-1° mean latitudinal shift towards the Equator and a mean increase in elevation of ~630m driven by 
an increase in maximum elevation. 95% of the permutation importance of the model was 
contributed to by minimum precipitation (82.5%), minimum temperature (5.9%), human influence 
index (3.0%), precipitation seasonality (2.5%) and surface roughness index (1.7%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A.F. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 





MaxTemp 





. — 3.0% 




-- 








MinTemi) 





"^N^ 5.9°/ 




— 








etpsLim 





1.4% 




— 















f 2.5% 




0 229 

ndvi 




^—^ 0-3°/ 


• 











,^<)%^^ 










._ 






0 1237 

RealMAR 




^ 0.2% 




._ 


— 






4j'r1 




radiation 






0.1% 










25.528 399 


Ci 13 


















0.4% 



E 



82.5% 



TempSea 




whpos 











0% 


- 



#62 - Moupin pika (Ochotona thibetana) 
n = 95 

Expert: Deyan Ge, Institute of Zoology, Chinese Academy 
of Sciences 

Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 2km/year (Similar ecology to O.roylei) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.93 


Omission rate 


0.13 


Sensitivity 


0.87 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.52 


True Skill Statistic 


0.86 



Summary: The Moupin pika's bioclimatic envelope is predicted to increase by 10% with a -2° 
mean latitudinal polewards shift and a mean increase in elevation of ~370m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
minimum precipitation (43.5%), maximum temperature (38.5%), temperature seasonality (7.7%), 
annual evapotranspiration (3.7%) and mean annual precipitation (3.1%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A.F. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#63 - Thomas's pika (Ochotona thomasi) 
n = 16 

Expert: Andrew Smith, Arizona State University 
Expert evaluation: Good 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 1 km/year (Similar ecology to O.koslowi) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.58 


True Skill Statistic 


0.99 



Summary: The Thomas's pika's bioclimatic envelope is predicted to decrease by 70% with a 
-1.5° mean latitudinal polewards shift and a mean increase in elevation of ~590m driven by an 
increase in maximum and minimum elevation. 95% of the permutation importance of the model 
was contributed to by minimum precipitation (56.6%), mean annual temperature (23.4%), minimum 
temperature (14.2%) and evapotranspiration (2.8%). 



a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A.F. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 






#64 - Turuchan pika (Ochotona turuchanensis) 
n = 30 

Expert: Andrey Lissovsky, Zoological Museum of Moscow 

State University 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 15km/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.93 


Omission rate 


0.13 


Sensitivity 


0.87 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.50 


True Skill Statistic 


0.86 



Summary: The Turuchan pika's bioclimatic envelope is predicted to increase by 50% with a -2° 
mean latitudinal polewards shift and a mean increase in elevation of ~70m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
minimum temperature (21.6%), minimum precipitation (19.9%), human influence index (15.2%), 
annual evapotranspiration (10.9%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Johnston, C.H. 
2008) 




#65 - European rabbit (Oryctolagus cunciculus) Model evaluation metric 

n = 22,712 

Expert: Neil Reid, Queen's University Belfast 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: lkm/year (Expert) 
Status: UNMODELLABLE 



AUC 


0.62 


Omission rate 


0.16 


Sensitivity 


0.84 


Specificity 


0.39 


Proportion correct 


0.62 


Kappa 


0.23 


True Skill Statistic 


0.23 



Summary: The European rabbit's bioclimatic envelope is predicted to increase by 30% with a -2° 
mean latitudinal polewards shift and a mean decrease in elevation of ~10m. 95% of the permutation 
importance of the model was contributed to by annual evapotranspiration (21.1%), mean annual 
temperature (19.0%), minimum temperature (13.9%), annual water balance (12.3%), mean annual 
precipitation (11.2%), human influence index (9.8%), precipitation seasonality (7.6%) and 
normalised difference vegetation index (2.9%). 
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#66 - Amami rabbit {Pentalagus furnessi) 
n = 9 

Expert: Fumio Yamada, Forestry and Forest Products 

Research Institute, Japan 

Expert evaluation: Good 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


0.99 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.95 


True Skill Statistic 


0.99 



Summary: The Amami rabbit's bioclimatic envelope is predicted to increase by 150% with a -1° 
mean latitudinal shift towards the Equator and a mean decrease in elevation of ~25m. 95% of the 
permutation importance of the model was contributed to by minimum temperature (97.3%), 
temperature seasonality (1.7%) and annual evapotranspiration (0.8%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Yamada, F & Sugimura, K. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#67 - Bunyoro rabbit (Poelagus marjorita) 
n = 8 

Expert: David Happold, Australian National University 
Expert evaluation: Poor 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 2km/year (Similar ecology to Pronolagus sp.) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.89 


True Skill Statistic 


0.99 



Summary: The Bunyoro rabbit's bioclimatic envelope is predicted to decrease by 90% with a -1° 
mean latitudinal shift towards the Equator and a mean increase in elevation of ~200m driven by an 
increase in minimum elevation.. 95% of the permutation importance of the model was contributed 
to by annual evapotranspiration (61.7%), temperature seasonality (31.4%) and number of months 
with a positive water balance (6.3%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Collins, K. & Smith, A.T. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#68 - Greater red rock hare {Pronolagus crassicaudatus) 
n = 7 

Expert: Kai Collins, University of Pretoria 
Expert evaluation: Poor 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 2km/year (Similar ecology to P.randensis) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.98 


Proportion correct 


0.98 


Kappa 


0.06 


True Skill Statistic 


0.98 



Summary: The Greater red rock hare's bioclimatic envelope is predicted to decrease by 65% with 
a -3° mean latitudinal polewards shift and a mean decrease in elevation of ~20m driven by a 
decrease in maximum elevation. 95% of the permutation importance of the model was contributed 
to by human influence index (31.7%), temperature seasonality (29.3%) solar radiation (16.9%) and 
minimum precipitation (18.3%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 




d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Matthee, C, et al. 2008) 
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#69 - Jameson's red rock hare (Pronolagus randensis) 
n = 27 

Expert: Kai Collins, University of Pretoria 
Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 
Dispersal distance: 2km/year (Expert) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.98 


Omission rate 


0.04 


Sensitivity 


0.96 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.55 


True Skill Statistic 


0.96 



Summary: The Jameson's red rock hare's bioclimatic envelope is predicted to decrease by 70% 
with a -5° mean latitudinal polewards shift and a mean increase in elevation of ~325m driven by an 
increase in maximum and minimum elevation. 95% of the permutation importance of the model 
was contributed to by mean annual temperature (46.2%), maximum temperature (23.3%), minimum 
precipitation (16.3%), minimum temperature (7.0%) and temperature seasonality (1.4%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Matthee, C, et al. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#70 - Smith's red rock hare (Pronolagus rupestris) 
n = 9 

Expert: Kai Collins, University of Pretoria 
Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 2km/year (Similar ecology to P.randensis) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.98 


Proportion correct 


0.98 


Kappa 


0.07 


True Skill Statistic 


0.98 



Summary: The Smith's red rock hare's bioclimatic envelope is predicted to decrease by 60% with 
a -5° mean latitudinal polewards shift and a mean decrease in elevation of ~60m. 95% of the 
permutation importance of the model was contributed to by temperature seasonality (44.0%), 
minimum precipitation (18.5%), mean annual precipitation (18.3%), normalised difference 
vegetation index (6.3%), human influence index (5.2%) and minimum temperature (3.4%). 




d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#71 - Hewitt's red rock hare (Pronolagus saundersiae) 
n = 9 

Expert: Kai Collins, University of Pretoria 
Expert evaluation: Poor 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 2km/year (Similar ecology to P.randensis) 
Status: UNMODELLABLE 

Summary: The Hewitt's red rock hare's bioclimatic envelope is predicted to decrease by 100% 
with a -1° mean latitudinal polewards shift and a mean increase in elevation of ~15m driven by an 
increase in minimum elevation. 95% of the permutation importance of the model was contributed to 
by temperature seasonality (28.7%), precipitation seasonality (24.3%), maximum temperature 
(18.0%), annual evapotranspiration (16.0%), human influence index (5.6%), minimum temperature 
(1.9%) and minimum precipitation (1.6%). 



Model evaluation metric 



AUC 1.00 

Omission rate 0.00 

Sensitivity 1.00 

Specificity 1.00 

Proportion correct 1 .00 

Kappa 1.00 

True Skill Statistic 1.00 




#72 - Volcano rabbit (Romerolagus diazi) 
n = 31 

Expert: Jose Antonio Martinez-Garcia, Universidad 
Autonoma Metropolitana, Mexico 
Expert evaluation: Poor 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Average for island species) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.79 


True Skill Statistic 


0.90 



Summary: The Volcano rabbit's bioclimatic envelope is predicted to decrease by 100% with a 
-0.2° mean latitudinal polewards shift and a mean increase in elevation of ~ 1500m driven by 
increases in minimum and maximum elevation. 95% of the permutation importance of the model 
was contributed to by temperature seasonality (75.2%), precipitation seasonality (16.2%) and 
minimum temperature (6.8%). 



a) ~ 1930s (1900-1949) 




.„ t 



b) ~ 1980s (1950-2000) 





c) ~ 2020s (2010-2039) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Malpica, F.J.,et al. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#73 - Swamp rabbit (Sylvilagus aquaticus) 
n = 66 

Expert: Robert Kissell, Memphis State University 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 25km/year (Expert) 

Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.09 


Sensitivity 


0.91 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.76 


True Skill Statistic 


0.91 



Summary: The Swamp rabbit's bioclimatic envelope is predicted to increase by 200% with a -2° 
mean latitudinal polewards shift and a mean increase in elevation of ~60m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
temperature seasonality (38.1%), mean annual temperature (36.7%), mean annual precipitation 
(12.7%), precipitation seasonality (4.7%) and minimum precipitation (3.6%). 




a) ~ 1930s (1900-1949) 




b) ~ 1980s (1950-2000) 




c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A.F. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#74 — Desert cottontail (Sylvilagus audubonii) 
n = 1040 

Expert: Consuelo Lorenzo, Departamento Conservacion de 
la Biodiversidad, Chiapas 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 7.5km/year (Similar ecology to S.palustris) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.94 


Omission rate 


0.08 


Sensitivity 


0.92 


Specificity 


0.96 


Proportion correct 


0.96 


Kappa 


0.78 


True Skill Statistic 


0.88 



Summary: The Desert cottontail's bioclimatic envelope is predicted to increase by 5% with a -1° 
mean latitudinal polewards shift and a mean increase in elevation of ~30m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
precipitation seasonality (36.6%), annual evapotranspiration (27.1%), mean annual temperature 
(13.6%), minimum temperature (9.1%), maximum temperature (3.0%), annual water balance 
(2.1%), minimum precipitation (1.9%) and human influence index (1.7%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#75 - Brush rabbit (Sylvilagus bachmani) 
n = 263 

Expert: Consuelo Lorenzo, Departamento Conservacion de 
la Biodiversidad, Chiapas 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Similar ecology to S.transitionalis) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.96 


Omission rate 


0.08 


Sensitivity 


0.92 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.89 


True Skill Statistic 


0.91 



Summary: The Brush rabbit's bioclimatic envelope is predicted to decrease by 15% with a -3° 
mean latitudinal polewards shift and a mean decrease in elevation of ~10m. 95% of the permutation 
importance of the model was contributed to by mean annual temperature (31.5%), minimum 
precipitation (24.5%), precipitation seasonality (14.4%), mean annual precipitation (7.1%), annual 
water balance (6.4%), minimum temperature (5.1%), normalised difference vegetation index (4.1%) 
and temperature seasonality (2.7%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J.et al. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#76 - Forest rabbit (Sylvilagus brasiliensis) 
n = 181 

Expert: Jorge Salazar-Bravo, Texas Tech University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 7.5km/year (Similar ecology to S.palustris) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.73 


True Skill Statistic 


0.89 



Summary: The Forest rabbit's bioclimatic envelope is predicted to decrease by 50% with a -6° 
mean latitudinal polewards shift and a mean increase in elevation of ~210m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by annual 
evapotranspiration (71.4%), temperature seasonality (11.6%), normalised difference vegetation 
index (4.9%) and minimum temperature (4.6%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 
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#77 - Manzano mountain cottontail (Sylvilagus cognatus) 

n = l 

Expert: Jennifer Frey, New Mexico State University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Similar ecology to R.diazi) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.52 


True Skill Statistic 


0.99 



Summary: The Manzano mountain cottontail's bioclimatic envelope is predicted to decrease by 
90% with a -2° mean latitudinal polewards shift and a mean increase in elevation of ~230m driven 
by an increase in minimum elevation. 95% of the permutation importance of the model was 
contributed to by annual water balance (40.8%), minimum temperature (21.7%), precipitation 
seasonality (11.3%), mean annual temperature (8.2%), temperature seasonality (6.1%), minimum 
precipitation (4.3%) and mean annual precipitation (2.6%). 




a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A. F. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 




#78 - Mexican cottontail (Sylvilagus cunicularius) 
n = 76 

Expert: Jorge Vazquez, Laboratorio de Ecologfa del 
Comportamiento, UAT-UNAM 
Expert evaluation: Medium 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 7.5km/year (Similar ecology to S.palustris) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.95 


Omission rate 


0.10 


Sensitivity 


0.90 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.73 


True Skill Statistic 


0.89 



Summary: The Mexican cottontail's bioclimatic envelope is predicted to decrease by 15% with a 
-0.5° mean latitudinal polewards shift and a mean increase in elevation of ~200m driven by an 
increase in minimum elevation. 95% of the permutation importance of the model was contributed to 
by temperature seasonality (48.3%), precipitation seasonality (42.8%), normalised difference 
vegetation index (3.5%) and minimum temperature (2.2%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#79 - Dice's cottontail (Sylvilagus dicei) 
n = 8 

Expert: Jan Schipper, Arizona State University 
Expert evaluation: Poor 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Similar ecology to S.cognatus) 
Status: UNMODELLABLE 



Model evaluation metric 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.73 


True Skill Statistic 


0.99 



Summary: The Dice's cottontail's bioclimatic envelope is predicted to decrease by 50% with a 
-1° mean latitudinal shift towards the Equator and a mean decrease in elevation of ~50m driven by a 
decrease in maximum elevation. 95% of the permutation importance of the model was contributed 
to by temperature seasonality (93.7%) and annual water balance (2.6%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A.F. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#80 - Eastern cottontail (Sylvilagus floridanus) 
n = 1104 

Expert: Jorge Vazquez, Laboratorio de Ecologfa del 
Comportamiento, UAT-UNAM 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 7.5km/year (Similar ecology to S.palustris) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


0.92 


Omission rate 


0.09 


Sensitivity 


0.91 


Specificity 


0.93 


Proportion correct 


0.93 


Kappa 


0.69 


True Skill Statistic 


0.84 



Summary: The Eastern cottontail's bioclimatic envelope is predicted to increase by 20% with a 
-2° mean latitudinal polewards shift and a mean increase in elevation of ~90m driven by an increase 
in minimum and maximum elevation. 95% of the permutation importance of the model was 
contributed to by annual evapotranspiration (29.5%), precipitation seasonality (26.2%), minimum 
temperature (11.5%), temperature seasonality (8.5%), minimum precipitation (8.1%), maximum 
temperature (3.3%), normalised difference vegetation index (2.9%), maximum precipitation (2.8%) 
and annual water balance (2.3%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#81 - Tres Marias cottontail (Sylvilagus graysoni) 
n = 6 

Expert: Consuelo Lorenzo, Departamento Conservacion de 
la Biodiversidad, Chiapas 
Expert evaluation: Good 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Average for island species) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


1.00 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


1.00 


Proportion correct 


1.00 


Kappa 


1.00 


True Skill Statistic 


1.00 



Summary: The Tres Marias cottontail's bioclimatic envelope is predicted to decrease by 20% 
with a no latitudinal polewards shift and a mean increase in elevation of ~25m driven by an increase 
in minimum elevation. 95% of the permutation importance of the model was contributed to by 
human influence index (56.0%), minimum precipitation (32.4%), precipitation seasonality (5.2%) 
and annual evapotranspiration (3.2%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#82 - Omilteme cottontail (Sylvilagus insonus) 

n = 3 

Expert: Alejandro Velazquez, UNAM-Canada 
Expert evaluation: Good 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Similar ecology to S.dicei) 
Status: MODELLABLE 



Model evaluation metric 



AUC 


1.00 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


1.00 


Proportion correct 


1.00 


Kappa 


1.00 


True Skill Statistic 


1.00 



Summary: The Omilteme cottontail's bioclimatic envelope is predicted to decrease by 80% with a 
no latitudinal polewards shift and a mean increase in elevation of ~120m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by 
normalised difference vegetation index (78.5%), precipitation seasonality (11.3%) and temperature 
seasonality (5.3%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Romero Malpica, F.J. & Rangel 
Cordero, H. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#83 - San Jose brush rabbit (Sylvilagus mansuetus) 
n = 9 

Expert: Tamara Rioja Pardela, Universidad de Ciencias y 
Artes de Chiapas, Mexico 
Expert evaluation: Good 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Average for island species) 
Status: MODELLABLE 



AUC 


1.00 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


1.00 


Proportion correct 


1.00 


Kappa 


1.00 


True Skill Statistic 


1.00 



Summary: The San Jose brush rabbit's bioclimatic envelope is predicted to decrease by 25% with 
a no latitudinal polewards shift and a mean increase in elevation of ~30m driven by an increase in 
minimum elevation. 95% of the permutation importance of the model was contributed to by 
temperature seasonality (36.7%), human influence index (34.0%), minimum precipitation (21.8%) 
and precipitation seasonality (3.2%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Lorenzo, C. & Alvarez- 
Castafieda, S. 2011) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#84 - Mountain cottontail (Sylvilagus nuttallii) 
n = 290 

Expert: Jennifer Frey, New Mexico State University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic and habitat 

Dispersal distance: 7.5km/year (Similar ecology to S.palustris) 
Status: MODELLABLE 



AUC 


0.95 


Omission rate 


0.09 


Sensitivity 


0.91 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.78 


True Skill Statistic 


0.90 



Summary: The Mountain cottontail's bioclimatic envelope is predicted to decrease by 20% with a 
~ 1 ° mean latitudinal polewards shift and a mean increase in elevation of ~40m driven by an increase 
in minimum elevation. 95% of the permutation importance of the model was contributed to by mean 
annual temperature (64.0%), maximum temperature (26.9%), temperature seasonality (2.8%) and 
minimum precipitation (1.4%). 






a) -1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Smith, A.T. & Boyer, A.F. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#85 - Appalachian cottontail (Sylvilagus obscurus) 
n = 39 

Expert: Michael Barbour, Auburn University 
Expert evaluation: Medium 
Data: Modern and historic 
Envelope: Climatic only 

Dispersal distance: O.Olkm/year (Similar ecology to S.dicei) 
Status: MODELLABLE 



AUC 


0.97 


Omission rate 


0.05 


Sensitivity 


0.95 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.73 


True Skill Statistic 


0.95 



Summary: The Appalachian cottontail's bioclimatic envelope is predicted to decrease by 10% 
with a -3° mean latitudinal polewards shift and a mean increase in elevation of ~70m. 95% of the 
permutation importance of the model was contributed to by annual evapotranspiration (42.0%), 
minimum temperature (24.2%), minimum precipitation (18.8%) and precipitation seasonality 
(10.0%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Barry, R. & Lazell, J. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#86 - Marsh rabbit (Sylvilagus palustris) 
n = 25 

Expert: Bob McCleery, University of Florida 
Expert evaluation: Good 
Data: Only modern 
Envelope: Climatic and habitat 
Dispersal distance: 7.5km/year (Expert) 
Status: MODELLABLE 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.75 


True Skill Statistic 


0.99 



Summary: The Marsh rabbit's bioclimatic envelope is predicted to increase by 90% with a -1° 
mean latitudinal polewards shift and a mean increase in elevation of ~2m driven by an increase in 
maximum elevation. 95% of the permutation importance of the model was contributed to by surface 
roughness index (73.1%), minimum precipitation (18.3%) and precipitation seasonality (4.9%). 






a) ~ 1930s (1900-1949) 



b) ~ 1980s (1950-2000) 



c) ~ 2020s (2010-2039) 





Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Faulhaber, C.A. & Smith, A.T. 
2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#87 - Robust cottontail (Sylvilagus robustus) 
n = 9 

Expert: Dana Lee, Oklahoma State University 

Expert evaluation: Poor 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: O.Olkm/year (Similar ecology 

Status: UNMODELLABLE 



AUC 


0.94 


Omission rate 


0.11 


Sensitivity 


0.89 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.27 


True Skill Statistic 


0.88 



S.dicei) 



Summary: The Robust cottontail's bioclimatic envelope is predicted to decrease by 90% with a 
-4° mean latitudinal polewards shift and a mean increase in elevation of ~480m driven by an 
increase in minimum elevation. 95% of the permutation importance of the model was contributed to 
by precipitation seasonality (68.5%), minimum precipitation (8.2%), temperature seasonality 
(7.1%), minimum temperature (4.5%), mean annual precipitation (3.6%), annual water balance 
(2.6%) and number of months with a positive water balance (1.6%). 






Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Ruedas, L. & Smith, A.T. 2008) 



d) ~ 2050s (2040-2069) 



e) ~ 2080s (2070-2099) 
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#88 - New England cottontail (Sylvilagus transitionalis) 
n = 18 

Expert: John Litvaitis, University of New Hampshire 

Expert evaluation: Medium 

Data: Modern and historic 

Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Expert) 

Status: MODELLABLE 



AUC 0.99 

Omission rate 0.00 

Sensitivity 1.00 

Specificity 0.99 

Proportion correct 0.99 

Kappa 0.68 

True Skill Statistic 0.99 



Summary: The New England cottontail's bioclimatic envelope is predicted to increase by 110% 
with a -1° mean latitudinal polewards shift and a mean increase in elevation of ~70m driven by an 
increase in maximum elevation. 95% of the permutation importance of the model was contributed 
to by annual evapotranspiration (23.5%), temperature seasonality (22.2%), minimum temperature 
(17.4%), normalised difference vegetation index (16.0%), mean annual temperature (8.6%), 
maximum temperature (5.9%) and minimum precipitation (5.9%). 




a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 




Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Barry, R., et al. 2008) 



d) ~ 2050s (2040-2069) e) ~ 2080s (2070-2099) 




#89 - Venezuelan lowland rabbit (Sylvilagus varynaensis) 
n = 6 

Expert: Daniel Lew, Venezuelan Institute of Scientific 
Research, Ecology Centre, Biodiversity Unit 
Expert evaluation: Poor 
Data: Only modern 
Envelope: Climatic and habitat 

Dispersal distance: 3km/year (Similar ecology to S.transitionalis) 
Status: UNMODELLABLE 



AUC 


0.99 


Omission rate 


0.00 


Sensitivity 


1.00 


Specificity 


0.99 


Proportion correct 


0.99 


Kappa 


0.92 


True Skill Statistic 


0.99 



Summary: The Venezuelan lowland rabbit's bioclimatic envelope is predicted to decrease by 
100% with a -1.5° mean latitudinal polewards shift and a mean increase in elevation of ~275m 
driven by an increase in minimum elevation. 95% of the permutation importance of the model was 
contributed to by temperature seasonality (97.7%). 



a) ~ 1930s ( 1 900- 1 949) b) ~ 1980s ( 1 950-2000) c) ~ 2020s (20 1 0-2039) 



Black = predicted suitable 
envelope 

Red polygon = IUCN range 
(Durant, P. & Guevara, M.A. 
2008) 






